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ABSTRACT 

We formulate a new, revised and coherent understanding of the structure and 
dynamics of the Large Magellanic Cloud (LMC), and its orbit around and interaction 
with the Milky Way. Much of our understanding of these issues hinges on studies of the 
LMC line-of-sight kinematics. The observed velocity field includes contributions from 
the LMC rotation curve V(R'), the LMC transverse velocity vector vt, and the rate of 
inclination change di/dt. All previous studies have assumed di/dt = 0. We show that 
this is incorrect, and that combined with uncertainties in vt this has led to incorrect 
estimates of many important structural parameters of the LMC. We derive general 
expressions for the velocity field which we fit to kinematical data for 1041 carbon stars. 
We calculate vt by compiling and improving LMC proper motion measurements from 
the literature, and we show that for known v t all other model parameters are uniquely 
determined by the data. The position angle of the line of nodes is = 129.9° ± 6.0°, 
consistent with the value determined geometrically by van der Marel & Cioni (2001). 
The rate of inclination change is di/dt = —0.37 ± 0.22 masyr -1 = —103 ± 61 
degrees/Gyr. This is similar in magnitude to predictions from iV-body simulations by 
Weinberg (2000), which predict LMC disk precession and nutation due to Milky Way 
tidal torques. The LMC rotation curve V(R') has amplitude 49.8 ± 15.9 kins" 1 . This 
is 40% lower than what has previously (and incorrectly) been inferred from studies of 
HI, carbon stars, and other tracers. The line-of-sight velocity dispersion has an average 
value a = 20.2 ± 0.5 kms -1 , with little variation as function of radius. The dynamical 
center of the carbon stars is consistent with the center of the bar and the center of the 
outer isophotes, but it is offset by 1.2° ± 0.6° from the kinematical center of the HI. 
The enclosed mass inside the last data point is Mlmc(8-9 kpc) = (8.7 ± 4.3) x 10 9 M , 
more than half of which is due to a dark halo. The LMC has a considerable vertical 
thickness; its V/a = 2.9 ± 0.9 is less than the value for the Milky Way's thick disk 
(V/a w 3.9). Simple arguments for models stratified on spheroids indicate that the 
(out-of-plane) axial ratio could be ~ 0.3 or larger. Isothermal disk models for the 
observed velocity dispersion profile confirm the finding of Alves &; Nelson (2000) that 
the scale height must increase with radius. A substantial thickness for the LMC disk is 
consistent with the simulations of Weinberg, which predict LMC disk thickening due to 
Milky Way tidal forces. These affect LMC structure even inside the LMC tidal radius, 
which we calculate to be r< = 15.0 ± 4.5 kpc (i.e., 17.1° ± 5.1°). The new insights into 
LMC structure need not significantly alter existing predictions for the LMC self-lensing 
optical depth, which to lowest order depends only on a. The compiled proper motion 
data imply an LMC transverse velocity vt = 406 kms -1 in the direction of position 
angle 78.7° (with errors of ~ 40kms~ 1 in each coordinate). This can be combined with 
the observed systemic velocity, v S y S = 262.2 ±3.4 kms" 1 , to calculate the LMC velocity 
in the Galactocentric rest frame. This yields wlmc = 293 ± 39 kms -1 , with radial 
and tangential components fLMC.rad = 84 ± 7 kms -1 and t>LMC,tan = 281 ± 41 kms -1 , 
respectively. This is consistent with the range of velocities that has been predicted by 
models for the Magellanic Stream. The implied orbit of the LMC has an apocenter 
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to pericenter distance ratio ~ 2.5 : 1, a perigalactic distance ~ 45 kpc, and a present 
orbital period around the Milky Way ~ 1.5 Gyr. The constraint that the LMC is 
bound to the Milky Way provides a robust limit on the minimum mass and extent of 
the Milky Way dark halo: M MW > 4.3 x 10 11 M and r h > 39 kpc (68.3% confidence). 
Finally, we present predictions for the LMC proper motion velocity field, and we 
discuss how measurements of this may lead to kinematical distance estimates of the 
LMC. 

Subject headings: galaxies: distances and redshifts — galaxies: kinematics and 
dynamics — Local Group — Magellanic Clouds. 
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1. Introduction 

The Large and Small Magellanic Clouds (LMC and SMC) are close companion galaxies of our 
Milky Way. Their interaction with the Milky Way has produced the Magellanic Stream, which 
spans more than 100° across the sky (e.g., Westerlund 1997; Wakker 2002). It consists of gas that 
trails the Magellanic Clouds as they orbit the Milky Way. A less prominent leading gas component 
was recently discovered as well (Lu et al. 1998; Putman 1998). Models of the Magellanic Stream 
hold the promise of providing important constraints on the mass and gravitational potential of 
the Milky Way dark halo (e.g., Lin, Jones & Klemola 1995). However, this requires an accurate 
understanding of the three-dimensional velocity of the Magellanic Clouds. While the systemic 
line-of-sight velocities of both the LMC and SMC have been accurately determined using Doppler 
shifts, a determination of their transverse velocities in the plane of the sky has proven to be much 
more difficult. We restrict our attention in the present paper to the LMC, which has been studied 
more and is understood better than the SMC. Several proper motion studies of the LMC have been 
published (Jones, Klemola & Lin 1994; Kroupa, Roser & Bastian 1994; Kroupa & Bastian 1997; 
Drake et al. 2002; Pedreros, Anguita & Maza 2002), but the uncertainties in the measurements are 
considerable. Accurate knowledge of the LMC transverse velocity has therefore remained elusive. 

The uncertainty in the transverse velocity of the LMC has not only restricted models of the 
Magellanic Stream, but also the interpretation of the observed LMC line-of-sight velocity field. 
As one moves away from the LMC center, the transverse velocity component of the LMC center 
of mass (CM) ceases to be perpendicular to the line of sight. This causes a spurious solid-body 
rotation component in the observed line-of-sight velocity field (Feast, Thackeray & Wesselink 
1961). Initially this effect was used to obtain a rough estimate of the transverse velocity of the 
LMC under the assumption (appropriate for a circular disk) that its kinematic line of nodes 
(defined as the line of maximum velocity gradient) should be equal to the photometric major axis 
(Feitzinger, Schmidt-Kaler & Isserstedt 1977; Meatheringham et al. 1988; Gould 2000). However, 
this assumption is now known to be invalid. Robust geometric methods show that the line of nodes 
of the LMC (defined as the intersection of the galaxy-plane and the sky-plane) does not coincide 
with the photometric major axis (van der Marel & Cioni 2001, hereafter Paper I), and hence, 
that the LMC is intrinsically elongated (van der Marel 2001, hereafter Paper II). The most recent 
analyses of the LMC velocity field have corrected for the spurious solid-body rotation component 
using the available LMC proper motion measurements (e.g., Kim et al. 1998; Alves & Nelson 
2000). However, the corrections involved are large (because the spurious component is comparable 
to the intrinsic LMC rotation component) and highly uncertain (because the LMC proper motion 
is not known accurately). In addition, all previous studies have assumed that the LMC viewing 
angles have no explicit time dependence, contrary to the results of A-body simulations which 
predict precession and nutation of the LMC disk due to Milky Way tidal torques (Weinberg 2000). 
All this probably explains why analyses of the LMC velocity field have yielded somewhat puzzling 
results. For example, Alves & Nelson (2000) obtained from a study of carbon-star velocities that 
the position angle of the kinematic line of nodes rises from 0kin = 143° ± 7° at a mean distance of 
4.5° from the galaxy center to 6 max = 183° ± 8° at 9.2°. The results obtained from HI kinematics 
by Kim et al. (1998) were qualitatively not very different. A kinematic twist of 40° ± 11° is quite 
unexpected, given that, over the same range of distances, the LMC has no twist in the stellar 



- 5 - 



number density contours (to an accuracy of ~ 3°; Paper II). More importantly, the inferred 
kinematic line of nodes differs considerably from the true line of nodes, for which a position angle 
= 122.5 ± 8.3° was determined geometrically from DENIS and 2MASS near-IR stellar catalogs 
(Paper I). While the LMC is not intrinsically circular, the implied difference between O and @kin 
is too large to plausibly attribute to non-circular orbits (Paper II). 

In the present paper we address these issues, and their consequences, by presenting the 
most detailed and sophisticated analysis to date of the LMC line-of-sight velocity field. In 
Section ^ we derive general expressions for the velocity fields of a rotating disk galaxy that 
has a non-negligible angular extent on the sky, both for the line-of-sight velocity component 
and for the velocity components in the plane of the sky. In Section ^ we present a conceptual 
discussion of the line-of-sight velocity field. We address which quantities are uniquely constrained 
by the observations, and which ones are not. In Section || we fit the general equations for the 
line-of-sight velocity field to the velocities available for 1041 carbon stars, and we discuss the 
results. In Section |H| we review what is known about the proper motion and distance of the LMC 
from other sources. In Section || we combine this knowledge with the information obtained from 
the analysis of the line-of-sight velocity field to determine di/dt, the rate at which the LMC 
inclination changes with time. In Section fj] we discuss the kinematical properties of the LMC disk 
as implied by our analysis. We address the position of the dynamical center, the disk rotation 
curve and the velocity dispersion profile, the position angle of the line of nodes, and the influence 
of non-circular orbits on the results of the analysis. In Section || we discuss the implications of the 
results for our understanding of the structure of the LMC, including its mass, tidal radius, scale 
height, dark halo and self-lensing optical depth. In Section |9| we discuss the implications of the 
results for the transverse velocity of the LMC, the three-dimensional velocity of the LMC in the 
Galactocentric rest frame, the orbit of the LMC around the Milky Way, the mass and extent of 
the Milky Way, and the origin of the Magellanic Stream. In Section [l(] we present predictions for 
the proper motion velocity field of the LMC, which may be observationally accessible with future 



astrometric missions. In Section 11 we discuss a new method for kinematical determination of the 
LMC distance, and its prospects for yielding results that are competitive with other methods. 
Conclusions are presented in Section 12. Appendix |A| discusses an improved analysis of the 



published proper measurements of Jones et al. (1994) and Pedreros et al. (2002), which were 
obtained for fields that are offset from the LMC CM. 



2. General Expressions for the Velocity Field of a Rotating Disk 

Many previous studies have given equations for the description of the velocity field of an 
external galaxy. However, it is customary to make the assumption that the angular size subtended 
by the galaxy is negligible. This is equivalent to assuming that the galaxy is at infinite distance 
and that 'the sky is flat' over the area of the galaxy. However, the LMC subtends a very large 
angle on the sky (more than 20° end-to-end) and this makes the usual equations inadequate for 
describing its kinematics. We therefore start by deriving general expressions for the velocity field 
of a rotating disk which are valid for a galaxy of arbitrary angular size. We address not only 
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the line-of-sight velocities, but also the velocity components perpendicular to the line of sight. 
Readers interested mostly in the application to the LMC may wish to skip directly to Section |3|. 



2.1. Coordinate Systems and Velocity Components 



We adopt the same coordinate systems as in Paper I. The position of any point in space is 
uniquely determined by its right ascension and declination on the sky, (a, 5), and its distance D. 
The point O with coordinates (ao,5o,Do) is the origin of the analysis, and is chosen to be the 
galaxy center of mass (CM). Angular coordinates (p,(ft) are defined on the celestial sphere, where 
p is the angular distance between the points (a, 5) and (ao> <5o)> and (ft is the position angle of the 
point (a, 5) with respect to (ao,<5o); see Figure |l[ In particular, (ft is the angle at (ao,5o) between 
the tangent to the great circle on the celestial sphere through (a, 5) and (ao,<5o), and the circle of 
constant declination 5o. By convention, (ft is measured counterclockwise starting from the axis that 
runs in the direction of decreasing right ascension at constant declination 5q. Equations (l)-(3) of 
Paper I allow (p,(ft) to be calculated for any (a, 5). 



At any given position (D, p, 
orthogonal components: 



a velocity vector can be decomposed into a sum of three 



Vl = 



dD 

~dt' 



v 2 



v 3 = Dsrnp—. 

at 



(1) 



Here, v\ is the line-of-sight velocity, and v 2 and V3 are the velocity components in the plane of 
the sky; see Figure ||. The goal of the present analysis is to obtain expressions for (vi,V2,v 3 ) as a 
function of (p, (j)). 

We introduce a Cartesian coordinate system (x,y,z) that has its origin at O, with the x-axis 
anti-parallel to the RA axis, the y-axis parallel to the declination axis, and the z-axis towards the 
observer; see Figures p] and |2|. The transformation equations from (D, p, cf>) to (x, y, z) are given 
by equation (5) of Paper I. The inverse transformation equations are: 



D 

P 



[x 2 + y 2 + (D 
arctan[(x 2 + 
arctan(y/x). 



z) 2 ] 1/2 , 

arctan[(x 2 + y 2 )/(D - z) 2 } 1/2 , 



(2) 



Substitution of equation (Q) in equation (|T|) yields after some manipulations 
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where (v x ,v y ,v z ) is the three-dimensional velocity in the (x,y,z) coordinate system. To describe 
the internal kinematics of the galaxy it is useful to adopt a second Cartesian coordinate system 
(a/, y', z') that is obtained from the system (x, y, z) by counterclockwise rotation around the z-axis 
by an angle 9, followed by a clockwise rotation around the new x'-axis by an angle i; see Figure |3|. 
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With this definition, the (x', y') plane is inclined with respect to the sky by the angle i (with 
face-on viewing corresponding to i = 0). The angle 9 is the position angle of the line of nodes 
(the intersection of the (x' , y')-plane and the (x,y)-plane of the sky), measured counterclockwise 
from the x-axis. In practice, i and 9 will be chosen such that the (x', y') plane coincides with the 
equatorial plane of the galaxy. The transformation equations from (x,y,z) to (x',y',z') are given 
by equation (6) of Paper I. The inverse transformations are 

1 sin i \ I x 1 \ 

(4) 

Upon taking the time derivative on both sides this yields transformation equations from (y' x ,v', v' z ) 
to (v x ,v y ,v z ). These can be substituted in equation (|3|), which after some further manipulations 
yields 

sinp cos(4> — 9) [sin pcos i s\n(<j> — 9) + cos psini] [sin p sin i sin(<p — 9) — cos pcos i] \ / v' x \ 
cospcos(0 — 9) [cos pcos i sin(0 — 9) — sin psin i] [cos p sin i sin(0 — 9) + sin pcos i] I vL I . (5) 
— sin(if> — 9) cosicos(0 — 9) sini cas(<j> — 9) I \ v' z J 

For the analysis of the line-of-sight velocity field one needs only v\, but to interpret proper 
motion measurements one needs to know in addition how the velocities V2 and i>3 relate to the 
local directions of West and North. The proper motions in these directions are defined as 

fi\V = — cos 5(da/dt), fj,^ = (d6/dt). (6) 

To obtain expressions for these proper motions we start with equations (l)-(3) of Paper I, which 
relate (p, <fi) to the right ascension and declination (a, S) on the sky. Upon taking the time 
derivative of these equations one obtains relations between dp/dt and d<p/dt on the one hand, and 
da/dt and dS/dt on the other hand. These can be solved to obtain: 

(l/D)x(- S[n J ~ COS l) ( V2 ). (7) 





The angle T determines the rotation angle of the (^2,^3) frame on the skyQ; see Figure |T[ It is 
determined by 

cosT = [sin 5 cos 5o cos (a — ao) — cos 5 sin 5o] /sinp, 

sinT = [cos So sin(a — ao)] / sinp. (8) 

In practice we restrict the discussion to a planar system with all tracers in the (x',y') plane. In 
this case the distance D to a tracer is related to the CM distance D$ through equation (8) of 
Paper I. Substitution of that equation yields 

[cos i cos p — sin i sinp sin(^ — 9)] . M l ■ 




Dq cos i 





lr The angle T is close to the position angle on the sky, $ = (f> — 90° (see Section below), but is not identical to 
it. The two are explicitly related to each other through the relation sinT = (cos So/ cos 8) sin"!?. 
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In practice one needs to express this in observable units, which is achieved by using the relation 
i/Do = {Oi/fkms -1 ]) / (4.7403885 D / kpc)} mas/yr. (10) 



2.2. Velocity Field Model 

For a planar system, the velocity of a tracer can be written as a sum of three components: 





CM 




(11) 



The first component, ucMj is the velocity contributed by the space motion of the CM. The second 
component, v pn , is the velocity contributed by time- variations in the viewing angles i and 8 of the 
disk plane. Such variations arise if the orientation of the disk plane is not static in an inertial 
frame. For the case of the LMC, this occurs naturally as a result of the precession and nutation 
of the disk plane induced by external tidal torques (Weinberg 2000). The third component, V{ n t, 
is the velocity contributed by the internal motion of the tracer in the disk plane. We proceed by 
deriving expressions for each of these contributions. 



2.2.1. Center- of- Mass Motion 



The velocity of the center of mass can be described by quantities v sys , vt and 6t, such that 
v sys is the systemic velocity of the CM along the line of sight (positive when receding), vt is the 
transverse velocity of the CM, and $t is the direction of the transverse motion on the sky; see 
Figure |||. This implies 
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From equations (||) and (12) one obtains for the observable velocity components 
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2.2.2. Precession and Nutation 



To assess the effect of changes in % and 6, consider a point that has fixed coordinates in the 
(x',y',z') frame. By taking the time derivative of equation (Q) one obtains: 
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For a planar disk (V = 0), one has from equation (2) of Paper II: 

/ x' \ „ . / cos z cos(</> — 0) 
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cos i cos p — sin i sin p sin(0 — 6 1 ) 
These equations can be substituted in equation (||), which yields after some manipulations 

Do sin p 



(di I dt) sin(</> — 9) [cos i cos p — sin i sin p sin(0 — 8)] 
(di I dt) sin(0 — 9) [— cos i sin p — sin i cos p sin((/> — 0)] 
(di/dt) sin(</> - 6»)[- sin i cos(0 - 6)] + (de/dt) [cos i] 



(14) 



(15) 



(16) 



cos 2 cos p — sin % sm p sm(0 — t 

i 

For future reference it is useful to note that the line-of-sight component of this vector simplifies to 

v i,pn = A) (di/dt) sin/) sin(<p — 9). (17) 
Note that the quantity d6/dt does not affect the line-of-sight velocity; it does affect the predicted 



proper motions, as discussed in Section 10 below. 



2.2.3. Internal Rotation 



We consider the case in which the mean streaming (i.e., the rotation) in the disk plane can be 
approximated as being circular. This is for simplicity only; there is no reason why the streamlines 



couldn't be non-circular, and why this couldn't potentially be important. Section discusses the 
extent to which this may influence our results for the LMC. For circular motion one has 



( sV{R')y'/R' 
-sV(R')x'/R! 
V o 



(18) 



where R' is the polar radius in the (x',y') plane, V(R') is the mean streaming velocity at radius 
R', and s = ±1 is the 'spin sign' that determines in which of the two possible directions the disk 
rotates. With help of equation (|l5|) one obtains 



R! 



Dq sin p [cos 2 i cos 2 



+ sin 2 ( ( />-fl)] 1 / 2 



cos i cos p — sin i sin p sin(0 — 9) 



(19) 
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and 



Substitution in equation (| 



cos 2 i cos 2 ( 
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(20) 



[cos 2 i cos 2 ((f> -9) + sin 2 (</> - 6)] 1 / 2 



— sin i cos(0 — 9) [cos i cos p — sin i sin p sin(</> — 9)] 
sin i cos(</> — 9) [cos i sin p + sin i cos p s'm(<j> — 9)] 
- [cos 2 i cos 2 (0 - 0) + sin 2 (<£ — 0)] 



(21) 



2.3. Position Angles 

In practice it is useful to employ positions angles $>, and Ot that are measured from North 
over East in the usual astronomical convention: 

$ = 90°, 9 = 0-90°, 9 t = 0t-9O°. (22) 

These are the angles that we use henceforth; see Figures ffl, [3| and The equations for the velocity 
field depend on the position angles (ft, 9 and 0t only through the differences eft — 8 and eft — 9t, for 
which 

^ - e = $ - e, ((> - 9 t = $ - e t . (23) 



3. Information Content and Degeneracies of the Line-of-Sight Velocity Field 

The model described in Section [T^ provides a closed expression for the line-of-sight velocity 
field of a planar disk with circular streamlines. Combination of equations (11, 13, 17, 21, 23) 
yields 

Ulos(P) = Wsyscosp + u t sinpcos($ - 6 t ) + D (di/dt) sinjOsin($ - 9) - s V r (i?')/sinicos($ - 0), (24) 
with 

. . cos z cos p - sinzsin / osin(<I> - 6) 

R=D smp f, f = - — — — — . 2 — T^vnT?- ( 25 ) 

[cos 2 zcos 2 (<l> - 6) + sm z ($ - 6)] 1 /2 

The quantities that feature in these equations have been defined previously in Section [||, but it is 
useful to recap them briefly. The velocity v\ os = v\ is the component of the velocity along the line 
of sight. The quantities (p, <E>) identify the position on the sky: p is the angular distance from the 
CM, and $ is the position angle with respect to the CM (measured from North over East). The 
quantities (v sys ,vt,Qt) describe the velocity of the CM in an inertial frame in which the sun is at 
rest: t> sys is the systemic velocity along the line of sight, vt is the transverse velocity, and t is 
the position angle of the transverse velocity on the sky. The angles (i, O) describe the direction 
from which the plane of the galaxy is viewed: i is the inclination angle (i = for a face-on disk), 
and is the position angle of the line of nodes. The velocity V(R') is the rotation velocity at 
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cylindrical radius R' in the disk plane. Dq is the distance to the CM, and / is a geometrical factor. 
The quantity s = ±1 is the 'spin sign' that determines in which of the two possible directions the 
disk rotates. In the following we will set s = +1, with the understanding that V(R') can be both 
positive and negative. 



The first two terms on the right hand side of equation (24) describe the line-of-sight 



component of the CM velocity of the LMC. The third term describes the line-of-sight component 
of the velocities induced by the precession and nutation of the LMC disk plane. The last term 
describes the line-of-sight component of the LMC rotation (in the limit p — > 0, the last term 
reduces to the rotation field model used by, e.g., Alves & Nelson (2000) and many other authors). 
The first term in equation (24) differs from v sys only by an amount of order p 2 . For p = 10° one 



has for the LMC v sys cos p ~ v sys — 4 kms -1 . Hence, the large size of the LMC has only a minor 
influence on the systemic velocity component of the observed velocity field. The second term in 
equation ( f24[) is proportional to sinp, which corresponds to a solid-body rotation component. At 
p = 10°, this component has an amplitude between 50-100 kms -1 , depending on the exact value 
of vt- This exceeds the projected amplitude of the intrinsic velocity field of the LMC (the fourth 
term in eq. [Q), which is only ~ 30-40 kms -1 . The third term also corresponds to a solid-body 
rotation component. Its amplitude is comparable to that of the second term if di/dt ~ vt/Do, 
which is in the range 1-2 mas/yr, depending on the exact value of vt- Expressed in more relevant 
units, this corresponds to ~ 275-550 degrees/Gyr. This is larger than what one would expect on 
the basis of iV-body simulations, but only by a factor of a few (Weinberg 2000). The third term 
in equation (p4j ) therefore should not be neglected, although this has in fact been done by all 
previous authors. 

To gain insight into the interplay between the different components that contribute to the 
observed velocity field it is useful to decompose the transverse velocity vector of the CM into a 
sum of two orthogonal vectors, one along the position angle of the line of nodes 0, and one along 
+ 90° ; see Figure |2[ These vectors have lengths 

v tc = v t cos(0 4 - 0), v ts = v t sin(0 t - 0), (26) 

respectively. The proper motion vector (pwi^n) of the CM (see eq. ||) can be similarly 
decomposed into a sum of vectors along and + 90°, respectively. These vectors have lengths 

p c = — pw sin + pn cos 0, p s = — pw cos — pjy sin 0. (27) 
The transverse velocity and proper motion components of the CM are related through 

v x = D Q p w , v y = D p N , v tc = D pc, v ts = Dqp s , (28) 
where v x and v y are the components of the transverse velocity in the (x, y, z) coordinate system 



defined in Section 2.1. In the following it is useful to define 

w ts = v ts + D (di/dt) = D [p s + (di/dt)] (29) 

and 

*los = V\ os - v sys cos p. (30) 
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The equation (24) for the line-of-sight velocity field then reduces to: 

^los = [wts ship] sin(<3? — 0) + [vt c sinp — f V(R') sini] cos(<I> — 0). (31) 

Along the line of nodes one has that sin(<I> — 0) = and cos(<3? — 0) = ±1. This yields 

vios (along) = ±[vtc sin p — V{Dq tanp) sini cos p]. (32) 

Perpendicular to the line of nodes one has that cos(3> — 0) = and sin(3> — 0) = ±1, and therefore 

^los (perpendicular) = ±wt s sin p. (33) 

This implies that perpendicular to the line of nodes v\ os is linearly proportional to sin p. By 
contrast, along the line of nodes this is true only if V(R') is a linear function of R' . This is not 
expected to be the case, because galaxies do not generally have solid-body rotation curves; disk 
galaxies tend to have flat rotation curves, at least outside the very center. This implies that, at 
least in principle, both the position angle of the line of nodes and the quantity wt s are uniquely 
determined by the observed velocity field: is the angle along which the observed v\ os are best 
fit by a linear proportionality with sin p, and Wt s is the proportionality constant. The systemic 
velocity, i> sys , and the position of the CM on the sky are also uniquely determined, because they 
represent points of symmetry in the velocity field. However, V(R'), i and vt c are not uniquely 
determined. Let V(R'), i and TJt c be the true values of these quantities. Then any combination of 
V(R'), i and vt c with 

V(R') = ( ^] V(R') + frfc-^) tan " (34) 
\ sini J sini 

for all R' will provide exactly the same predictions both along and perpendicular to the line of 
nodes. While the predictions may be very subtly different at intermediate angles, the quality 
of realistic datasets will be insufficient to discriminate between different models that satisfy 
equation (j34|). Hence, the component Vf C of the transverse velocity is unconstrained by the 
observed velocity field, unless we assume some prior knowledge about the shape of the rotation 
curve V(R'). Conversely, the rotation curve V(R') cannot be determined unless vt c is known. 



4. The LMC Line-of-Sight Velocity Field Traced by Carbon Stars 

The kinematical properties of the LMC have been studied using many different tracers, 
including HI (e.g., Kim et al. 1998), star clusters (Freeman, Illingworth & Oemler 1983; Schommer 
et al. 1992), planetary nebulae (Meatheringham et al. 1988), HII regions and supergiants 
(Feitzinger et al. 1977), and carbon stars. The latter have yielded the largest and most useful 
datasets in recent years, and we therefore restrict our analysis to carbon star data. We included 
and merged two different radial velocity data sets. The first is the one obtained by Kunkel et 
al. (1997), which was analyzed previously by Alves & Nelson (2000). This data set covers mostly 
the periphery of the LMC, and is distributed fairly homogeneously in position angle. We excluded 
all (inter-cloud) stars further than p = 13° from the LMC center, because it is not clear to what 
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extent they belong to the LMC and participate in its kinematics. This leaves a total of 468 stars. 
The second data set is the one obtained by Hardy, Schommer & Suntzeff (2002, in preparation), 
which was analyzed previously by Graff et al. (2000). The 573 stars in this data set were selected 
from Blanco et al. (1980), and their spatial distribution reflects the distribution of fields in that 
survey. The stars sample the inner parts of the LMC, but with a discontinuous distribution in 
radius and position angle. The combined dataset of 1041 stars samples the entire area of the LMC, 
but not homogeneously. The same data set, although with subtly different selection criteria, was 
studied previously by Alves (2000); the spatial distribution in the LMC of these carbon stars with 
known line-of-sight velocities is shown in his figure 1. 



We fitted the velocity field given by equation ( pip to the data, by minimizing the RMS 
residual of the fit using a downhill-simplex routine (Press et al. 1992). The position of the CM 
and the parameters (v sys , vt s , Vt c , ©) were all treated as free parameters that were optimized in the 
fit. The inclination i is essentially degenerate with the amplitude of the rotation curve V(R'), 
cf. equation (^). It was therefore, without loss of generality, fixed to the value 

i = 34.7° ± 6.2° (35) 

determined in Paper I. The rotation curve V(R') was parameterized as 

v W' v °WTw (36> 

This corresponds to a velocity that increase as a power-law until some scale radius i?o, after which 
it flattens to the constant value Vo- The parameters (Vo,Ro,r]) were also optimized in the fit. The 
adoption of a parameterized form for V(R') artificially removes some of the degeneracy described 
by equation (^); while infinitely many functions V(R') can fit the data equally well, only some of 
these can be expected to be properly described by the adopted parameterization. The implications 
of this are addressed below. 

Once the best-fitting model has been identified, we calculate error bars on the model 
parameters using Monte-Carlo simulations. Many different pseudo data sets are created that are 
analyzed similarly as the real data set. The dispersions in the inferred model parameters are a 
measure of the formal 1-a error bars on the best-fit model parameters. Each pseudo data set is 
created by calculating for each star a velocity v\ os that is the sum of a random Gaussian deviate 
Av and the velocity predicted by the best-fit model. The Gaussian dispersion of the Av is chosen 
to correspond to the RMS residual of the best-fit model. 

The data and the best-fitting model are shown in Figure ||. For all radii 1° < p < 9° the fit is 
about as good as could be expected. At small radii, p < 1°, there appears to be some systematic 
discrepancy between the data and the model. However, the data in the central degree (top left 
panel of Figure [5]) provide too sparse coverage of the full range of position angles to draw any 
conclusions from this. There also appear to be minor discrepancies in the data-model comparison 
at radii p > 9°. This could be due to tidal distortions in the LMC disk and its velocity field at 
large radii. Either way, these discrepancies do not affect the primary results of this paper in any 
significant way; fits to restricted radial ranges yielded identical results. 
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The right ascension a and declination 5 of the CM and the parameters v sys , wt s and G are all 
determined unambiguously. The RMS residual of the fit is the (globally averaged) line-of-sight 
velocity dispersion a of the LMC. The model fit yields: 

a CM = 5 h 27.6 m ± 3.9 m , 

S CM = -69.87° ±0.41°, 

v sys = 262.2 ±3.4 km s~\ 

Wts = -402.9 ± 13.0 kms -1 , 

G = 129.9° ±6.0°, 

a = 20.2 ±0.5 km s -1 , (37) 

where the CM position is given in J2000.0 coordinates. 

As discussed in Section H the quantity vt c and the rotation curve V(R') are not uniquely 
determined by the line-of-sight velocity field. The parameterized model that best fits the data 
has vtc = 599 kms -1 and the rotation curve V(R') that is shown as a solid curve in Figure ^. 
However, every other vt c will provide an equally acceptable fit, provided that V(R') is modified 
according to equation (|34|). The alternative rotation curves thus obtained are shown as dashed 
curves in Figure ^ for several values of vt c - The value of vt c is constrained by the line-of-sight 
velocity field only through the fact that some of the rotation curves V(R') needed to fit the data 
are unplausible. If vt c is high, then V(R') is implausibly large for a low-luminosity galaxy such 
as the LMC. The LMC is slightly fainter than M33, My = -18.5 vs. My = -18.9, respectively 
(van den Bergh 2000), and the LMC and M33 have similar disk scale lengths (Paper II; Regan & 
Vogel 1994). According to the Tully-Fisher relation, the circular velocity of the LMC should be 
lower than for M33. For M33, V(R') m 120 kms" 1 at R' = 11.6 kpc (Corbelli & Salucci 2000), the 
largest radius included in our study of the LMC. Figure || shows that for Vf c = 600 kms -1 , the 
LMC rotation curve levels off at V(R') ~ 130 kms -1 . This makes models with vt c ^ 600 kms -1 
quite implausible. On the other hand, if vt c is low, less than ~ 200 kms -1 , then V(R') has a sign 
change within the radial range out to which its starlight can be traced. This is never observed in 
disk galaxies. Based on these considerations we conclude that 

v tc 6 [200,600] kms -1 . (38) 

The results described in this section were found to be very robust. We experimented with 
fits to subsets of the data, fits over restricted radial ranges, fits with different parameterizations 
for V(R'), and different algorithms for the determination of the parameters that characterize the 
velocity field. All of this yielded results that are statistically equivalent to those described above. 



5. The Proper Motion and Distance of the LMC 

The analysis of the LMC line-of-sight velocity field has provided constraints on the velocities 
vtc and wt s (eqs. [37| and |38| ). These velocities are directly related to the proper motion vector 
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{nwy MJv) and distance Do of the LMC CM through equations (§7]), ( |28| ) and (|29|). To make further 
progress it is therefore useful to review what is already know about these quantities. 



5.1. Proper Motion 

Table [l] lists the proper motion measurements that are available for the LMC. They are 
from the following sources: Kroupa et al. (1994), using stars from the PPM Catalogue; Jones 
et al. (1994), using photographic plates with a 14 year epoch span; Kroupa & Bastian (1997), 
using Hipparcos data; Drake et al. (2002), using data from the MACHO project^]; and Pedreros 
et al. (2002), using CCD frames with an 11 year epoch span. The measurements of Jones et 
al. (1994) and Pedreros et al. (2002) pertain to fields in the outer parts in the LMC disk. These 
measurements require corrections for the orientation and rotation of the LMC disk. The proper 
motions listed in the table are based on the corrections that we derive in Appendix [A], which are 
more accurate than the corrections derived in the original papers. The table does not include 
the measurement published by Anguita, Loyola & Pedreros (2000). Their result for the proper 
motion in the declination direction (jUtv = 2.9 ± 0.2 masyr -1 ) is highly (10a) inconsistent with all 
other studies. This probably indicates that the measurement contains some unidentified error; see 
Pedreros et al. (2002) and Drake et al. (2002) for some discussion of this issue. 

Figure shows the data from Table |l| in the (/j-w^n) plane. The ellipses around the data 
points indicate the corresponding 68.3% confidence regions. The proper motion measurements are 
in fair agreement with each other, in the sense that they fall in the same part of the diagram. The 
weighted averag eg of the proper motion measurements in Table |] is 

m = -1.68 ±0.16 mas yr _1 , ^ N = 0.34 ± 0.16 masyr" 1 . (39) 

This is the LMC proper motion that we will use in the subsequent discussion. The corresponding 
68.3% confidence region is indicated in Figure |?] by a heavy ellipse. 



5.2. Distance 

The distance Dq to the LMC is a very important quantity in the calibration of the cosmological 
distance ladder. Most of the methods that are customarily used to study the Hubble constant of 



2 The MACHO proper motion measurement is based on a preliminary analysis of the data. However, the result of 
the final analysis is not expected to be considerably different (Drake 2002, priv. comm.) 

3 The 1-a error in a weighted average of N measurements Zi is generally given by E^=i l/(Az;) 2 ] -1//2 . However, 
this is correct only if all the errors in the data are Gaussian random errors. In the present case there is some evidence 
that this is an oversimplification. For example, the measurements by Jones et al. (1994) and Pedreros et al. (2002) 
are mutually inconsistent at the ~ 3a level. In taking the weighted average we therefore increased all the errors by 
25%. This causes the x 2 that measures the residuals of the data with respect to their weighted average to be equal 
to the number of degrees of freedom. 
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the Universe are ultimately calibrated by the distance to the LMC. The period-luminosity relation 
of Cepheids provides the most important, but by no means the only example. Of course, there now 
exist methods that provide information on the Hubble constant without requiring knowledge of the 
distance to the LMC. These use, e.g., gravitational lensing, the Sunyaev-Zeldovich effect, or the 
Cosmic Microwave Background radiation. However, the systematic uncertainties in these methods 
have not yet been reduced to levels where they can replace the more traditional techniques. The 
distance to the LMC therefore remains a topic of intense interest and debate in modern astronomy. 

Many techniques have been used to estimate the LMC distance. Unfortunately, there continue 
to be systematic differences between the results from different techniques that exceed the formal 
errors. It is beyond the scope of the present paper to review this topic in detail. Instead, we refer 
the reader to the recent reviews by, e.g., Westerlund (1997), Gibson et al. (2000) and Freedman et 
al. (2001). Here we follow Freedman et al. by adopting 

m-M = 18.50 ±0.10 (40) 

for the LMC distance modulus, where m — M = 51og(Do/ kpc) + 10. The implied distance is 
Dq = 50.1 ±2. 5 kpc. At this distance, a proper motion of lmasyr -1 corresponds to 238 ±12 km s -1 
(eq. [K|), and 1 degree on the sky corresponds to 0.875 ± 0.044 kpc. 

New and improved distance determinations for the LMC continue to become available at a 
rapid rate. Recent measurements that were not yet included in the compilation of Freedman et 
al. (2001) include, e.g., measurements using the Tip of the Red Giant Branch (Cioni et al. 2000), 
eclipsing binaries (Fitzpatrick et al. 2002), RR Lyrae variables (Benedict et al. 2002) and the Red 
Clump (Alves et al. 2002, in prep.). The results from these recent studies are all consistent with 
the value adopted in equation (pCf). 



6. Precession and Nutation of the LMC Disk 



Equation (29) can be solved for di/dt to obtain 

di/dt = (w ts /D ) - fi s . (41) 

The quantity wt s is determined by the line-of-sight velocity field, as is the line-of-nodes position 
angle G (eq.|37]]). The proper motion component fi s is determined by the literature average in 
equation (|39[), with help of the definition in equation (p7[). The distance is given by the literature 
average in equation (00). The rate of inclination change di/dt is therefore uniquely determined by 



the available data. Evaluation of equation (41) yields 

di/dt = -0.37 ± 0.22 masyr" 1 = -103 ± 61 degrees/ Gyr, (42) 

where a simple Monte-Carlo scheme was used for propagation of errors. 

The LMC is the first galaxy for which it has been possible to measure di/dt. It is therefore 
useful to ask whether the inferred value is plausible, given our understanding of the LMC and 
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its orbit around the Milky Way. Weinberg (2000) performed iV-body simulations to assess the 
influence of the Milky Way on the structure and dynamics of the LMC. He found that the tidal 
torques from the Milky Way are expected to induce precession and nutation in the symmetry 
axis of the LMC disk plane. His figures 3 and 4 show the variations in the directions of the 
LMC symmetry axis for a generic simulation of the LMC (i.e., not fine-tuned to reproduce all 
the currently observed features of the LMC). On average, the rate of directional change in the 
three-dimensional angle is ~ 60°/ Gyr (i.e., 0.22 masyr -1 ). However, there are considerable 
variations with time, and at certain times the rate is several times larger. Whether this directional 
change is seen by an observer as a change in inclination or as a change in the line-of-nodes position 
angle depends on the exact location of the observer. Either way, the observed rate di/dt given by 
equation (^) has the same order of magnitude as the generic predictions by Weinberg (2000). We 
therefore interpret the observed di/dt as an indication of precession and nutation induced by tidal 
torques from the Milky Way. 

The line-of-sight velocity field does not determine the second derivative d 2 i/dt 2 . It also does 
not provide any information on time variations in the position angle of the line of nodes. It is 
therefore not possible to calculate either the past or the future variation of the orientation of the 
LMC symmetry axis in an inertial frame tied to the Milky Way. If one assumes that d 2 i/dt 2 is 
always equal to zero, then the LMC would make one revolution every 3.5 ±2.1 Gyr. However, this 
number provides no real insight. The simulations of Weinberg show that LMC symmetry axis is 
expected to undergo nutation (i.e., oscillation with d 2 i/dt 2 ^ 0) and not tumbling motion (i.e., 
full revolutions). 



7. Kinematical Properties of the LMC disk 

7.1. Dynamical Center 

Previous authors who have modeled the dynamics of tracers in the LMC have fixed the 
dynamical center a priori, generally to coincide with the kinematical center of the HI gas. Kim 
et al. (1998) find this HI kinematical center to be at a = 5 h 17.6 m and 5 = -69°1' (J2000.0 
coordinates), with an accuracy that is probably better than ~ 0.2° in each coordinate. It has 
been a well-known result for many years that this position differs by almost a full degree from 
the center of the LMC bar (e.g., Westerlund 1997). The data from 2MASS yield for the latter 
a = 5 h 25. l m ± 0.1 m and 5 = -69° 47' ± T (Paper II). Our analysis in Section | is the first to 
provide an accurate measurement of the dynamical center of the stars in the LMC, by leaving it as 
a completely free parameter in the fit. We obtain «cm = 5 h 27.6 m ±3.9 m and 5cm = — 69°52'±25' 
(eq. |^7| ). This is consistent with the position of the center of the bar. It is also consistent with 
the position of the center of the outer isophotes of the LMC (corrected for the effect of viewing 
perspective), which is at approximately a = 5 h 29 m and 6 = — 69°30' (Paper II). However, the 
dynamical center of the stars differs from the dynamical center of the gas by as much as 1.2° ±0.6°. 
These results suggest that the gas kinematics are quite disturbed, while the stellar kinematics 
show little evidence for peculiar behavior. This is not particularly surprising in view of other 
knowledge about the gas in the Magellanic Clouds, which is quite disturbed in general (witness 
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the Magellanic Bridge and the Magellanic Stream). These results make it unlikely that any study 
of the HI in the LMC will ever lead to an accurate understanding of its structure and dynamics. 



7.2. Rotation Curve 

For studies of the dynamics and mass distribution of the LMC it is useful to have an 
unparameterized description of its rotation and velocity dispersion profiles. To this end we have 
performed fits of the velocity field given by equation (p4|) to individual rings in the plane of the 
LMC disk. The rotation velocity V(R') and the line of sight position angle O were assumed to be 
constant over each ring, and were determined so as to best fit the data. The resulting kinematical 
profiles are listed in Table |2| and are shown in Figure [f| In the fits we kept the quantities 
(aCM, <5cm, v sys , wts, i) fixed to their previously determined values (eqs. |35| and |37]]). The value of 
Vtc was fixed to that implied by the LMC proper motion; from equations (27, |2^, |37|, B9L ^) one 
obtains 

vtc = 253 ±52 km s" 1 . (43) 

This is not inconsistent with the weak constraint vt c £ [200,600] kms" 1 that was obtained from 
the analysis of the line-of-sight velocity field (eq. [38Q. 

The inferred rotation curve V(R') rises linearly in the central region and is roughly flat at 
V « 50 kms -1 for R 1 ^ 4 kpc (i.e., R'/D ^ 0.08). The error bars in V(R') do not take into 
account the errors in (acM> ^CM)f sys , Wt s , i, Vtc), which were all kept fixed. In reality, the errors 
in i and v% c do add additional uncertainty to the rotation curve. The uncertainty in i (eq. j35|) 
causes a ~ 15% uncertainty in the normalization of V(R'), cf. equation (p4|). No other part of our 
analysis depends on the exact choice of the inclination. The influence of the 52 kms -1 error in vt c 
(eq. |43|]) can be assessed visually from Figure ||. 

The innermost data point in Figure || at R'/Dq 0.01 (i.e., R' « 0.5 kpc) is clearly discrepant 
from the data points at larger radii, both in terms of its kinematic position angle and its 
rotation velocity V . Although this could be attributed to non-circular streaming motions in the 
region of the bar, this result has only very low significance. The top left panel of Figure || shows 
that the data in this region provide only very sparse coverage of the full range of position angles. 
Better data are needed to address the kinematics in the central kpc with any confidence. 

Alves & Nelson (2000) previously analyzed the carbon-star kinematics of the LMC. 
Like most authors, they didn't use the fully correct expression for the velocity field of 



a rotating disk at finite distance (last term in eq. [24Q and they assumed di/dt = 0. 
They corrected for the transverse motion of the LMC using an assumed proper motion 
(HWiH-n) = (—1.65 ± 0.20,-0.17 ± 0.22) masyr -1 . This is an average of the proper motion 
measurements by Jones et al. (1994) and Kroupa & Bastian (1997). This proper motion, combined 
with di/dt = 0, implies v tc = 325 ± 57 kms^ 1 and w ts = -220 ± 62 kms" 1 (cf. eqs.||, |8], |?], ||, 



40). The value of wt s is quite different from the value that we have obtained here from a fit to the 



velocity field, w% a = — 402. 9± 13.0km s 1 (cf. eq. (37j). The result of enforcing an incorrect value of 
wts on the solution is that one obtains an incorrect value for the position angle of the line of nodes. 
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Alves & Nelson obtained a kinematic line of nodes that is both twisting and inconsistent with the 
geometrically determined line of nodes (see Section |l|). Since the rotation curve is determined 
by the velocities measured along the line of nodes, one also obtains an incorrect estimate of the 
rotation curve. Alves & Nelson obtained a rotation curve amplitude of ~ 70 kms" 1 , which exceeds 
ours by ~ 40%. 

The analysis by Kim et al. (1998) of the HI velocity field suffered from similar shortcomings. 
They corrected for the transverse motion of the LMC using the proper motion advocated by Jones 
et al. (1994), (hw,^n) = (—1.37 db 0.28, —0.18 ±0.27) masyr" 1 . This proper motion, combined 
with di/dt = 0, implies v% c = 276 d= 69 kms -1 and wt s = —175 db 72 kms" 1 . Again, the value of 
Wf S is quite different from the value that we have obtained here from a fit to the velocity field. 
So, like Alves & Nelson, Kim et al. obtained a kinematic line of nodes that is both twisting and 
inconsistent with the geometrically determined line of nodes, as well as a rotation curve amplitude 
that exceed ours by ~ 40%. 



7.3. Velocity Dispersion Profile 

The velocity dispersion profile that we derive for the carbon stars in the LMC is shown in 
the middle panel of Figure ^. It is not very different from the velocity dispersion profile derived 
by Alves & Nelson (2000). The dispersion is 20-22 kms -1 between 1 and 3.5 kpc from the center, 
followed by a decline to 16-17 kms -1 between 3.5 and 7 kpc from the center, and a subsequent 
increase to 21-22 kms -1 between 7 and 9 kpc. 

The different stellar populations in the LMC are not characterized by a single velocity 
dispersion. As in the Milky Way, younger populations have a smaller velocity dispersion than 
older populations. A summary of measurements for various populations is given by Gyuk, Dalai 
k, Griest (2000). They range from a rj 6 kms" 1 for the youngest populations (e.g., supergiants, 
HII regions, HI gas) to a ~ 30 kms" 1 for the oldest populations (e.g., old long-period variables, 
old clusters). Any discussions in the remainder of the paper that use our velocity dispersion 



measurements, e.g., those in Section B.3 concerning the LMC scale height, apply strictly only to 
carbon stars. On the other hand, the carbon stars are part of the intermediate-age population 
which is believed to be fairly representative for the bulk of the mass in the LMC. In this sense, the 
results inferred for the carbon star population are believed to be generic for the LMC as a whole. 



7.4. Line of Nodes Position Angle 

The fit to the velocity field described in Section |3] implies a line-of-nodes position angle 
O = 129.9° dz 6.0° (eq. p7[j ). This agrees within the errors with the geometrically determined 
line-of-nodes position angle, = 122.5 dz 8.3° (Paper I). This agreement was not in any way built 
into the model, and it therefore provides an important consistency check on the validity of our 
approach. This is particularly important because it resolves the discrepancy that emerged from 
previous kinematical analyses of the LMC (see Section |l]). It also provides further confirmation 
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of the finding from Paper II that the position angle of the line nodes is very different from the 
position angle of the LMC major axis (PA ma j = 189.3° ± 1-4°), and hence that the LMC is not 
intrinsically circular. 

The bottom panel of Figure [] shows that the radial profile of is approximately flat (with 
the exception of the innermost point, which has very low significance; see Section [7^). The 
constancy of O as a function of radius should come as no surprise. Our velocity field model does 
not allow for radial twists in 0. The best-fitting value for wt s is therefore, by definition, the 
value that yields the least possible radial variation in O. One might wonder if this is not the 
equivalent of forcing a square peg into a round hole. After all, the projected morphology of the 
LMC shows a considerable isophotal twist between the region of the bar and the outer parts of 
the disk (Paper II). While this is certainly a worry, we have found little evidence that there is a 
corresponding twist in the kinematic line of nodes. First, the results of the present paper, which 
assume an absence of kinematic twists, provide a highly consistent view of the LMC. Second, there 
is no evidence for any isophotal twist at projected radii p > 4.5°, suggesting that at least at these 
radii there are probably no twists in the kinematic line of nodes. Kinematical fits restricted to the 
range p > 4.5° yielded results that were statistically consistent with those for the full radial range. 



7.5. Influence of Non-Circular Streamlines 

It was shown in Paper II that the LMC is not intrinsically circular, but instead has an 
ellipticity e = 1 — (b/a) = 0.31 ± 0.01. However, the gravitational potential of a mass distribution 
is always rounder than the mass distribution, which provides some a priori justification for the use 
of a model in which the streamlines are circular. There is also some a posteriori justification from 
the fact that the kinematically inferred position angle of the line of nodes, ©kin = 129.9° ± 6.0° 



(eq. [37|), agrees with the geometrically determined value, O = 122.5 ± 8.3° (Paper I). Detailed 
dynamical models of elliptical disks generally predict a (small) misalignment between ©km an d © 
(e.g., Schoenmakers, Franx & de Zeeuw 1997). The fact that no such misalignment is detectable 
with the available statistics (©kin — © = 7.4° ± 10.2°) suggests that a model with circular 
streamlines may not be unreasonable in the present context. Note, however, that this argument 
cannot be reversed. The fact that ©kin — © is small does not by itself imply that the streamlines 
must be nearly circular (Jalali & Abolghasemi 2002). 

Unfortunately, there are no unambiguous theoretical methods for calculating the velocity field 
of an elliptical disk. Some results have been obtained for gaseous systems (Jalali & Abolghasemi 
2002), but stellar systems are more complicated because they can have an anisotropic velocity 
dispersion (i.e., pressure) tensor. It is therefore difficult to estimate the systematic errors in our 
results due to the simplifying assumption of circular streamlines. It should be noted though, that 
even if the streamlines are not circular, this will not invalidate our approach at a basic level. Even 
in an elliptical disk there is generally a position angle T along which the line-of-sight velocity 
component is zero (Schoenmakers et al. 1997). Unlike in a circular disk, one will not generally 
have that T = + 90°, where is the true geometrical line of nodes. However, along the position 
angle T one still has that v\ os = ±wt s sin p, as in equation (^) . So for an elliptical disk one can 
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still, at least in principle, constrain T and wt s from the line-of-sight velocity field. However, the 
actual algorithm employed here will not be quite sufficient to yield highly accurate results. We 
rely on fitting the velocity field given by equation (24), which for an elliptical disk is not generally 
valid at all position angles <£. It is not clear how important this is for the present analysis, but it 
will certainly become more important as larger samples of carbon star line-of-sight velocities will 
become available in the future. This will decrease the formal error random errors in the modeling, 
at which point deviations from circular streamlines may start to play a more important role. 



8. Structure of the LMC Disk 

An important result of the present study is that we find the rotation curve amplitude of the 
LMC to be ~ 40% lower than suggested on the basis of previous analyses (see Section 7.2). This 
changes our understanding of several important structural parameters and properties of the LMC. 



8.1. Mass 

Figure |6| shows that the rotation curve of the LMC is approximately flat from R ~ 4 kpc out 
to the last data point at R = 8.9 kpc. The weighted mean of the rotation velocity data points 
in this radial range is V = 49.8 kms -1 . There are several uncertainties in this rotation curve 
amplitude: (i) an error of ±2.1 kms -1 due to the finite number of data points; (ii) an error 
of ±7.8 kms -1 due to the error in our knowledge of the LMC inclination (eq. |]35]|); and most 
importantly (iii) an error of approximately ±13.7 kms -1 due the error in our knowledge of vt c 
(eq. |43] and Figure ||). Addition of these errors in quadrature yields V = 49.8 ± 15.9 kms -1 . The 
weighted mean line-of-sight velocity dispersion over the same radial range is a = 16.9 ± 0.5 kms -1 . 

Due to asymmetric drift, the average azimuthal streaming velocity V is always smaller than 
the velocity V c \ rc of a tracer on a circular orbit in the equatorial plane. For the solar neighborhood 
in the Milky Way, observations indicate that V 2 iTC « V 2 + 5.5<7jj (Dehnen & Binney 1998). For 
the LMC we expect similarly that = V 2 + kct 2 , where a is the observed line-of-sight velocity 
dispersion and k is a model dependent factor. Since the formal error in V is quite large, there 
is no particular benefit in making very sophisticated models to determine k. A simple model 
that is adequate in the present context considers the equatorial plane of an axisymmetric system 
with an isotropic velocity distribution embedded in an isothermal dark halo. The equations of 
van der Marel (1991) then show that — k is equal to the logarithmic slope of the mass density in 
the equatorial plane. For an exponential disk this yields k = R/Rj, where Rd is the exponential 
disk scale length. From Paper II we know that the LMC is modestly well described by an 
exponential profile with Rj ~ 1.5 kpc. This implies that at the last measured data point n ~ 6 
and V circ = 64.8 ± 15.9 kms" 1 (the error on V c \ rc does not include a theoretical error on k, given 
that its influence would most likely be negligible compared to the formal error in V). 

The total mass of the LMC inside the last measured data point is Mlmc(-R) = ^Kh-c/^ 
where the gravitational constant G = 4.3007 x 10 -6 kpc( kms -1 ) 2 M© -1 . This yields 
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M LM c(8.9 kpc) = (8.7 ± 4.3) x 10 9 M . Despite its relatively large error bar, our estimate is 
considerably more accurate in a systematic sense than the results of previous studies. Nonetheless, 
we find that (after correction for the fact that different authors tend to quote masses enclosed 
within different radii) our result is similar to most other LMC mass estimates in the literature 
(e.g., Feitzinger et al. 1977; Meatheringham et al. 1988; Kim et al. 1998). The reason for 
this is that while previous authors generally inferred a higher rotation curve amplitude V (see 
Section [7.2j ), their general neglect of the asymmetric drift correction led to similar estimates of 
V^i rc . Our mass estimate does differ considerably from that presented by Schommer et al. (1992). 
They applied a statistical mass estimator to the velocities of a sample of 16 old LMC clusters to 
infer a mass Mlmc ~2x 10 10 M within an effective radius of ~ 6 kpc. This result is inconsistent 
with the analysis presented here. Several iV-body simulation studies of the LMC have assumed, 
following Schommer et al. , that ~ 2 x 10 10 Mq is enclosed inside the visible body of the LMC 
(e.g., Gardiner, Sawa & Fujimoto 1994; Gardiner & Noguchi 1996; Weinberg 2000). It is likely 
that these simulations have overestimated the self-gravity of the LMC compared to the Milky Way 
tidal force. 



8.2. Tidal Radius 

To calculate the tidal radius of the LMC we consider the gravitational accelerations at a point 
on the line that connects the Milky Way to the LMC. Let the point be at a distance (3Dq from 
the LMC center and at a distance [1 — (3]Dq from the Milky Way center. We assume that [3 <C 1; 
the calculations can easily be carried out to full accuracy, but this does not lead to significantly 
different results. Let omw be the gravitational acceleration due to the Milky Way and 3lmc the 
gravitational acceleration due to the LMC. The tidal radius corresponds to the point at which 
omw([1 - 0\Dq) - a M w(D ) = -a LM c((3D ). 

The traditional assumption in the calculation of the tidal radius is that both bodies can be 
approximated to be point masses. In this case one obtains that 

ftidai = Pi = [M LMC /2M MW ] 1/3 • (44) 

This equation was used, e.g., by Weinberg (2000; his appendix A). However, this equation is not 
really appropriate, because there is strong evidence that the mass of the Milky Way continues to 
rise linearly to beyond the distance of the LMC (e.g., Kochanek 1996; Wilkinson & Evans 1999). 
If one assumes instead that the Milky Way has a flat rotation curve with circular velocity Vo, 
while still assuming that the LMC is a point mass, then 

Aidai = (3 2 = [M LMC /Mmw(A))] 1/3 . (45) 

where Mmw(-Co) is the mass of the Milky Way inside a sphere of radius Dq. This yields a tidal 
radius that is a factor 2 1 / 3 = 1.26 larger than the value given by equation ([H|). If one assumes 
that the LMC itself can also not be approximated as a point mass, but instead has a constant 
circular velocity V c \ rc (i.e., mass rising linearly with radius), then 



Aidal = $3 = V c i TC /Vo- 



(46) 
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The kinematical data that are available for the LMC do not extend as far out as the 
tidal radius, so we do not know whether the LMC circular velocity curve at the tidal radius 
is flat or Keplerian. The LMC tidal radius, r t = AidaiAh must therefore be in the range 
#2 A) <; n <; foDo- To evaluate equation (ph we use M LM c = (8-7 ± 4.3) x 10 9 M from 



Section 8.1 and M MW (A)) = (4.9 ± 1.1) x 10 11 M from Kochanek (1996), which yields 



/?2 = 0.26 db 0.05. To evaluate equation (EH) we use V c \ rc = 64.8 ± 15.9 kms 1 from Section 8.1 and 



V = 206 ± 23 kms" 1 based on the mass given by Kochanek, which yields (3$ = 0.31 ± 0.08. We 
combine these results to obtain for the LMC tidal radius that rt = 15.0 ± 4.5 kpc. This corresponds 
to an angle on the sky of 17.1° ± 5.1°. 

The result thus calculated can be compared to our observational knowledge of the outer 
parts of the LMC. Data from the 2MASS survey show that the isopleths of RGB and AGB stars 
can be traced out to a major axis radius of ~ 9° and are regular out to this radius (Paper II). 
APM measurements of photographic Schmidt plates have revealed regular isopleths out to ~ 11.5° 
(Irwin 1991). Observations therefore indicate that rt 11.5°, consistent with value calculated 
above. Having said this, it must be kept in mind that the tidal radius is not a particularly 
well-defined quantity, especially not for a disk system such as the LMC. The Milky Way center 
does not lie in the plane of the LMC disk (Paper II). As a consequence, for points in the LMC 
disk the gravitational accelerations omw and 3lmc are not colinear. So there is no radius in the 
disk where they add up to zero. Also, as the LMC orbits the Milky Way, the angle between the 
acceleration vectors changes. The effect of the Milky Way tidal force therefore need not lead to a 
sharp truncation of the disk, but instead may lead to thickening of the disk and the formation of 
an extended halo of unbound particles (Weinberg 2000). 



8.3. Scale Height 

The observed velocity dispersion profile contains important information on the scale height 
of the LMC as function of radius, as discussed in detail by Alves & Nelson (2000). They found 
that the velocity dispersion profile does not fall steeply enough with radius to be consistent with a 
constant scale- height disk. From this they concluded that the disk must be 'flared', with a scale 
height that increases radially outward. This can arise naturally as a result of tidal forces from the 
Milky Way. These forces become relatively more important (compared to the LMC self-gravity) 
as one moves to larger radii. Indeed, the iV-body simulations by Weinberg (2000) predict a 
considerable thickness for the LMC disk as a result of the Milky Way tidal forces (although, 
somewhat surprisingly, the vertical velocity dispersion in his simulations remains small). 

It is straightforward to refit the isothermal flared disk model of Alves & Nelson (2000; their 
eq. [37]) to our new kinematical results and preferred LMC disk parameters. For an isothermal 
disk the vertical density profile is proportional to sech 2 (z/zo), where zq can vary with disk radius. 
The input to the model of Alves & Nelson (2000) includes the total disk mass (Section |S.lD , 
the radial exponential scale length of the disk (Paper II), the total mass of the LMC dark halo 
(Section ^J|), the mean density of the Milky Way's dark halo at the LMC distance (Section |9.4[) 
and approximate corrections for the finite thickness and finite radial extent of the disk. The best 
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fit to our velocity dispersion data (Figure |6| and Table ^) yields zq = 0.27 kpc at the LMC center, 
rising to zq = 1.5 kpc at a radius of 5.5 kpc. If one prefers to model the disk's thickness with an 
exponential instead of a sech 2 function, then the exponential scale height is approximately half of 
zq (Gilmore, King & van der Kruit 1990). 

It is useful to compare the LMC disk to the disk components of our Milky Way galaxy^. 
The thin disk of the Milky Way has V/a « 9.8 and the thick disk has V/a ps 3.9. For the LMC 



we have, averaged between 4 and 9 kpc as in Section 8.1, that V/a ~ 2.9 ± 0.9. The Milky Way 
thick disk has an exponential vertical scale height ~ 1 kpc (Binney & Merrifield 1998). Its 
exponential radial scale length is ~ 3 kpc, if one assumes that the Milky Way thin and thick 
disks have the same radial profile. Hence, at the solar radius (~ 3Rd), the Milky Way thick disk 
is characterized by a relative thickness Zd/Rd ~ 1/3. This is quantitatively similar to the values 
inferred for the LMC, upon application of Alves & Nelson's (2000) flared disk model to our new 
data. At 3Rd ~ 4.5 kpc, the model predicts zq = 1.1 kpc. The corresponding exponential scale 
height (see above) is Zd ~ 0.5 kpc so that Zd/Rd ~ 1/3. So both in terms of V/a as well as in 
terms of Zd/Rd, the LMC is very similar to the thick disk of the Milky Way and very different 
from the thin disk of the Milky Way. 

There are different approaches to estimate the actual shape and vertical scale height of the 
LMC disk. A useful alternative to the approach employed by Alves & Nelson (2000) is to assume 
that the LMC is spheroidal with axial ratio q (for reference, an infinitely thin disk has q <C 1 
whereas the very flattest elliptical galaxies have q ~ 0.3; Binney & Merrifield 1998). The axial 
ratio is directly related to the quantity V/a through hydrostatic equilibrium. For a self-gravitating 
system one can use the relation provided by the tensor virial theorem (Binney & Tremaine 1987). 
However, this relation is not accurate for the LMC, which is embedded in a dark halo (see 
Section below). So we consider another simple model instead: we assume that the LMC is 
embedded in a spherical isothermal halo, and that it has an isotropic velocity dispersion that is 
approximately constant as function of radius. The equations of van der Marel (1991) then yield 
that 

9 I 1 -1/2 

-(V/a) 2 ) , (47) 

K I 



2 + 



where, as in Section p.l| , — n is the logarithmic slope of the mass density in the equatorial plane: 
k ~ R/Rd with Rd ~ 1.5 kpc. This simple model suggest that the axial ratio ranges from q ~ 0.31 
at R = 3 kpc to q ~ 0.46 at R = 9 kpc. These results are qualitatively consistent with those from 
Alves & Nelson's (2000) flaring disk model. 

The LMC has V/a considerably larger than unity, and it therefore remains quite justified to 
regard the LMC as a rotationally supported disk-like system. However, it is most definitely not 
a thin disk. Although the axial ratio values that we have estimated are based on simple models, 
and need to be interpreted with care, it is clear that the LMC is quite thick. While this does 



4 We use the kinematical properties of the Milky Way disks as listed in Table 10.4 of Binney & Merrifield (1998). 
The velocity dispersion components o\r, and a z were used to calculate the velocity dispersion a that would be 
observed if the Milky Way were viewed with the same inclination as the LMC. 
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not appear to have been fully appreciated in the literature, it need not come as a surprise. The 
inferred thickness is entirely consistent with the predictions from simulations (Weinberg 2000). 



8.4. Dark Halo 

The total 5-band magnitude of the LMC is 0.63 (de Vaucouleurs & Freeman 1973), the 
B — V colorQ is 0.61 ± 0.03 (Bothun & Thompson 1988), and the V-band extinction is Ay « 0.4 
(Zaritsky 1999). This implies a total intrinsic luminosity Ly = 3.0 x 10 9 L©. Realistic stellar 
population synthesis and evolution models for disk galaxies predict a strong correlation between 
color and mass-to-light ratio. To fit the extinction-corrected B — V color, the mass-to-light ratio 
of the stellar population should be M/Ly « 0.9 ± 0.2 (Bell & de Jong 2001), with some residual 
dependence on the assumed initial mass function and its low-mass cut-off. This implies a total 
mass for the visible disk of the LMC of M ~ 2.7 x 10 9 M . The mass of the neutral gas in 
the LMC has been estimated to be 0.5 x 10 9 M© (Kim et al. 1998). The combined mass of the 
visible material in the LMC is therefore insufficient to explain the dynamically inferred mass 
M LM c(8.9 kpc) = (8.7 ± 4.3) x 10 9 M Q (Section f$A§. Consequently, the LMC must be embedded 
in a dark halo. This is consistent also with the fact that the observed rotation curve amplitude is 
relatively flat as a function of radius. To determine the properties of the LMC dark halo one must: 
(i) model the contributions to the circular velocity from gas, stars and dark matter; (ii) calculate 
the asymmetric drift; (iii) perform a detailed data-model comparison. Such a detailed analysis is 
beyond the scope of the present paper. Alves & Nelson (2000) did perform such modeling, but 
their results will need to revised in view of the present paper: they used an observed rotation 
curve amplitude that is 40% higher than the one that we have inferred here, and they ignored 
asymmetric drift. 

The observed rotation velocity V(R) listed in Table ^ rises very slowly; it does not reach its 
flat part until R ~ 4 kpc. This is different from the behavior of the circular velocity curves V c i TC (R) 
predicted by realistic models, which reach velocities near their maximum already at R ~ 2 kpc 
(Alves & Nelson 2000). This indicates, independent of model details, that V < V c i TC in the central 
region of the LMC. Our favored explanation for this is that it is the result of asymmetric drift. In 
the central R <^ 2 kpc, V/cr is of order unity. We do not think that it is due to some error in the 
rotation curve V{R) inferred from the observations. At small radii, V(R) is rather insensitive to 
errors in vt c , cf. Figure ||. There also is no strong dependence on the assumed position angle of 
the line of nodes. Alves (2000) also found that V < V c hc at R <^ 4 kpc, despite the use of a very 
different line of nodes position angle. Nonetheless, one must not forget that the central region 
of the LMC contains a strong bar. So it is quite likely that the streaming velocity field in this 
region is considerably more complicated than the circular streaming model that we have used in 
our analysis. A proper understanding of the kinematics and mass distribution of the LMC in its 
central few kpc will require detailed modeling of both bar dynamics and asymmetric drift. 



5 The discussion in Section 3c of Bothun & Thompson (1988) argues that the B — V colors of the LMC and SMC 
in their Table 1 are erroneously reversed. 
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8.5. Microlensing Optical Depth 

The observed microlensing optical depth towards the LMC is r Q bs = x 10 -8 with an 
additional 20-30% of systematic error (Alcock et al. 2000). To interpret this result it is critically 
important to know the self-lensing optical depth of the LMC, r se if. This subject was addressed 
and reviewed most recently by Gyuk et al. (2000). For their most favored set of LMC model 
parameters^ they predict r se if = 2.2 x 10 -8 , a factor of 5.5 less than the observed value. It is useful 
to address whether our new results on LMC structure change this conclusion. 

To lowest order, the LMC self-lensing optical depth depends exclusively on the observed 
velocity dispersion and not separately on either the galaxy mass or scale height (Gould 1995). 
The values that we have inferred here for the line-of-sight velocity dispersion of the LMC are 
not very different from those used by previous authors. This argues that, to lowest order, our 
new findings do not change the predictions for the LMC self-lensing optical depth. In principle, 
the expression derived by Gould (1995) is somewhat idealized. It relies on the assumptions 
that: (a) the LMC disk is thin; and (b) the observed velocity dispersion tracks the scale height 
as in simple equilibrium models. The first assumption is suspect based on our discussion in 
Section |8.3| and the second assumption has been shown to be violated in the iV-body simulations 
of Weinberg (2000). One could therefore argue that it is prudent to consider more sophisticated 
models. However, this was done by Gyuk et al. (2000), and they did not find results that differ 
greatly from the simple predictions obtained directly from the velocity dispersion. Alves Sz Nelson 
(2000) incorporated the flaring of the LMC disk, and this too was found not to increase r se if very 
significantly. The MACHO survey fields used to derive r b s are mostly situated in the central 2 
kpc of the LMC (Alcock et al. 2000), and at these radii the flaring is not yet substantial enough 
to make a big difference (see Section |8.3[ ). In view of these arguments, it appears unlikely that the 
revised understanding of LMC structure that has emerged from the present study will significantly 
increase the predicted LMC self-lensing optical depth. 

For LMC self-lensing models to become a viable explanation for the observed lensing events 
they will not only have to correctly predict r b s , but also the observed spatial and time-scale 
distributions of the lensing events (Gyuk et al. 2000; Alcock et al. 2000). Proposed tests of the 
lensing population based on the spatial distribution of the events may need to be revisited in light 
of our new and coherent understanding of LMC structure. For example, the tests discussed by 
Alcock et al. (2000; their section 5.3; see also Gyuk et al. 2000) are based on an incorrect line of 
nodes position angle of 170°. This is quite different from the actual value, which has now been 
found using several different methods to be in the range 120°-130° (Papers I and II; Section ||). 
This error could become important, particularly if larger samples of microlensing events become 
available. The flare of the LMC disk is also important for self-lensing models. It causes the spatial 
distribution to more resemble that expected for LMC halo lenses. The event rate will fall off more 
gradually with radius than if the disk had a constant thickness, and the near-to-far side event rate 
asymmetry will be less. It may be worth investigating new tests of the spatial event distribution 



6 The relevant estimate to use for comparison to the Alcock et al. (2000) data is the 30-field average in Table 4 of 
Gyuk et al. (2000). 



that build on our improved understanding of the LMC structure (e.g., a test based on the event 
distance perpendicular to the line-of-nodes). 



9. The Orbit of the LMC 

9.1. Systemic Velocity 

As mentioned in Section |4], the kinematical properties of the LMC have been studied 
previously using many different tracers. Each of these studies has yielded an estimate of the 
systemic velocity of the LMC. The result that has most often been adopted in models of the Milky 
Way dark halo and the Magellanic Stream (e.g., Wilkinson & Evans 1999; Gardiner et al. 1994) is 
the value v sys = 274 kms -1 inferred from HI data by Luks & Rohlfs (1992). This value is similar 
to the value f; sys = 279 kms -1 inferred subsequently, also from HI data, by Kim et al. (1998). 
However, the HI distribution and kinematics in the LMC are quite disturbed. For example, Luks 
& Rohlfs required two separate HI velocity components to fit their data, and as discussed in 
Section [O], the HI dynamical center probably does not coincide with the LMC CM. The true 
systemic velocity, i.e., the line-of-sight velocity of the CM, is therefore better estimated from 
discrete tracers. From our analysis of 1041 carbon stars[] we have obtained t> sys = 262.2 ± 3.4km s _1 
(eq. |37]]). The error in the determination of t> sys scales as l/V^/V, where N is the number of 
tracers. Previous studies with other tracers generally had much smaller values of N than available 
for the present study. For example, Meatheringham et al. (1988) studied 94 planetary nebulae; 
Schommer et al. (1992) studied 83 star clusters; etc. We conclude that our estimate of v sys is the 
most accurate one currently available. 



9.2. Transverse Velocity 

The observed proper motion of the LMC quoted in equation (|39|) yields an estimate of the 



transverse velocity of the LMC. With the help of equations (pq) and (40) one obtains 



Vx = -399 km s" 1 , v y = 80 km s" 1 , (48) 

where the x and y directions point towards the West and North, respectively. This result 
corresponds to a velocity vt = 406 kms -1 in the direction of position angle Qt = 78.7°. The 
uncertainty in the distance Dq causes the error Av in the transverse velocity to be larger in the 
direction Gj than in the direction Gj + 90°: 

Av(G t ) = 44kms-\ Av(O t + 90°) = 37 kms" 1 . (49) 



7 We verified that the two subsets of our carbon star data set, i.e., those obtained by Kunkel et al. (1997) and 
Hardy, Schommer & Suntzeff (2002, in preparation), respectively, yield mutually consistent values for v sys (to within 
the error bars). So there is no evidence to doubt the velocity calibration accuracy of either data set. 
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The 68.3% confidence region on the transverse velocity of the LMC is shown as a heavy ellipse in 
Figure || 

The analysis of the line-of-sight velocity field also constrains the transverse velocity of the 
LMC. This velocity is determined completely by the orthogonal components vt c and Vt s , measured 
along and perpendicular to the line of nodes, respectively (Figure ||). The component vt c is 
weakly constrained by the velocity field, vt c £ [200,600] kms -1 (eq. fl38| ). The component 
vts = Wis — Do(di/dt) (eq. [^]) can be calculated only if wt s , Dq and di/dt are all known. While wt s 
and Dq are rather tightly constrained, di/dt is unfortunately completely unknown. So in Section || 
we used this relation to estimate di/dt, using the value of vt s implied by the proper motion data. 
For the purpose of visualization it is now useful to reverse this situation for a moment, and assume 
that di/dt is known. In particular, let us assume that di/dt = 0. In this case one obtains from 
the observed line-of-sight velocity field an estimate of the transverse velocity of the LMC that is 
completely independent from the available proper motion and distance measurements. The 68.3% 
confidence region thus obtained is the trapezoid^ shown with dashed lines in Figure ||. It overlaps 
(albeit slightly) with the 68.3% confidence region obtained from the proper motion data. The 
agreement can be improved by varying di/dt, which shifts the trapezoid along the vt s direction 
shown in Figure ^. The agreement is optimized for di/dt = —0.37 mas yr™ 1 , the value derived in 
Section ||. The 68.3% confidence region for this value is shown as a solid trapezoid in Figure [8|. 
Note, that the shift in the (v x ,v y ) plane induced by this di/dt is fairly small. In addition, it was 
shown in Section || that the simulations of Weinberg (2000) actually predict a somewhat smaller 
absolute value of di/dt. So theory predicts that the true transverse velocity of the LMC should 
not be very different from the value predicted from the velocity field under the assumption that 
di/dt = 0. This is exactly what the proper motion measurements indicate. 



9.3. Three-Dimensional Space Motion 

Our knowledge of the transverse velocity of the LMC is becoming steadily more accurate as 
the number of independent proper motion measurements increases. We also know now, based 
on the results shown in Figure ||, that the proper motion measurements yield a picture that is 
consistent with the results from analysis of the line-of-sight velocity field. Therefore, there is no 
reason to believe that the proper motion measurements are plagued by large systematic errors. It 
is therefore prudent to ask what the measurements imply for the three-dimensional space motion 
of the LMC. This issue has been addressed previously by several authors, but the results from the 



Some brief explanation is in order as to why the confidence region on the transverse velocity obtained from the 
velocity field is a trapezoid. To calculate v x and v y it is insufficient to know the components v tc and v ta along and 
perpendicular to the line of nodes. One needs to know instead the components along fixed position angles on the sky. 
Let these angles be Go and Go + 90°, and let the corresponding velocity components be VtcO and VtsO- We choose 
Go = 129.9°, which is the best estimate of the line-of-nodes (eq. J37|). Propagation of errors then yields for the formal 
error in vt s o that Avmo = [(Avt s ) 2 + « 2 c0 (AG) 2 ] 1//2 . Here AG = 0.10 is expressed in radians and Av ta = 13.0 kms -1 
(eq. p^l). Because VtcO > 200kms _1 (eq. pq |) one finds that the error in vt s o is dominated by the error in the position 
angle of the line of nodes, and that Avtso ~ 0.10 VtcO- Visual representation of this relation yields a trapezoid. 
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present paper now allow us to address this with higher confidence. 

To determine the three-dimensional motion of the LMC with respect to the Milky- Way we 
need to correct for the reflex- motion of the sun. Following Gardiner et al. (1994) and Kroupa & 
Bastian (1997) we adopt a Cartesian coordinate system (X,Y,Z), with the origin at the Galactic 
Center, the Z-axis pointing towards the Galactic North Pole, the X-axis pointing in the direction 
from the sun to the Galactic Center, and the Y-axis pointing in the direction of the sun's Galactic 
Rotation. The position and velocity vectors of the sun are 

r = (-Ro, 0, 0), v Q = (U Q , V + V G , W Q ). (50) 

The quantities (i?o,Vo) are the distance of the sun from the Galactic Center, and the circular 
velocity of the Milky Way at the position of the sun. We use the standard IAU values Rq = 8.5 kpc 
and Vq = 220 kms" 1 (Kerr & Lynden-Bell 1986). The vector (Uq,V q ,W & ) is the sun's velocity 
with respect to the Local Standard of Rest for which we use the recent determination by Dehnen 
& Binney (1998): (U Q , V Q , W Q ) = (10.0 ± 0.4, 5.2 ± 0.6, 7.2 ± 0.4) kms" 1 . 

Let (l,b) be the Galactic coordinates of the LMC. The right ascension and declination of the 
dynamical center of the LMC determined in Section |4| yield, using equation (2.1) of Binney & 
Merrifield (1998)0, I = 280.531° and b = -32.523°. The unit vector from the sun in the direction 
towards the LMC is U\ = (cos b cos I, cos b sin /, sin b), so that the position vector of the LMC is 
fLMC = ^0 + DqU\. For the distance Dq to the LMC we adopt the value from equation (p0|). This 
implies 

f hMC = (-0.78, -41.55, -26.95) kpc. (51) 

By numerical differentiations of the unit vector u\ with respect to right ascension and declination, 
respectively, one can calculate the unit vectors ui and U3 in the directions West and North (in the 
plane of the sky, as seen from the sun). The unit vectors thus obtained, 

711 = (0.15410, -0.82897, -0.53764), 

u 2 = (0.06860, 0.55180, -0.83115), (52) 
u 3 = (-0.98567, -0.09120, -0.14190), 

form an orthonormal basis set of the (X, Y, Z) space. The velocity vector of the LMC corrected 
for the reflex motion of the sun is 



VLMC = Ve+ WsysWl + v x u 2 + v y u 3 , (53) 



where t> sys is the systemic velocity of the LMC (eq. |57j) and v x , v y are the transverse velocity 
components of the LMC (eqs. p^j49[ |). The total velocity of the LMC in the (X,Y,Z) system is 
vlmc = \vlmc\- The radial velocity component is t>LMC,rad = 4mc • rLMc/ln.Mcl- Equation 
yields for the distance of the LMC form the Galactic Center 

J-LMC = |r L Mc| = 49.53 kpc. (54) 



An error in this equation was corrected; the longitude of the North Celestial Pole is Icp = 122.932° . 



- 30 - 



The tangential velocity component of the LMC velocity is t>LMC,tan = { v lmc ~ ^LMCrad) 1 ^ 2 - These 
equations yield for the three-dimensional space velocity of the LMC: 

v LMC = (-56 ±36, -219 ±23, 186 ±35) kms" 1 , 

vlmc = 293 ±39 km s" 1 , t> LM c,rad = 84 ± 7 kms -1 , t>LMC,tan = 281 ± 41 kms" 1 . (55) 

The errors in equation (|55|) were obtained by propagation of the errors in v SJS , v x and v y using 
a simple Monte-Carlo scheme. The errors in the model parameters Rq, Vq, Uq, Vq, W & , Dq, 
and the position of the LMC center on the sky were not explicitly taken into account. However, 
it was verified that variation of these parameters within their error ranges does not change any 
of the values in equation (]55| ) by more than a few km/s. The only exception is Vq, which adds 
directly into the ^-component of vlmc- This changes fLMd ^LMC.rad and VLMC.tan as well. So a 
considerably non-standard value of Vo (e.g., Vq = 184 kms -1 , as advocated by Oiling & Merrifield 
1998) would change the results in equation (j55|) by a few tens of kms -1 . 



The values in equation (55) are not hugely dissimilar from those calculated by, e.g., Kroupa 
& Bastian (1997) and Pedreros et al. (2002). The main difference is that our error bars are 
considerably smaller. This is partly due to the use of a larger data set, and partly due to a more 
accurate method of error propagation. The results for the LMC in equation ( |55| ) now provide one 
of the highest accuracy space velocities of any Galactic satellite (velocities for other satellites can 
be found in, e.g., Wilkinson & Evans [1999] and Dinescu et al. [2001]). 



9.4. Mass and Extent of the Milky Way 

Models generally suggest that the Magellanic Stream represents material torn from the 
Magellanic Clouds during the previous perigalactic passage (e.g., Lin et al. 1995). The fact that 
this is not the first passage implies that the LMC must be bound to the Milky Way (but see Shuter 
1992 for an alternative view). So the binding energy £ = ^> — ^v 2 must be positive (Binney Sz 
Tremaine 1987). It is useful to study what this implies for the mass and extent of the Milky Way. 

We adopt a simple model in which the Milky Way dark halo has a spherically symmetric mass 
density p(r) with an isothermal profile and a sharp cut-off at r = r^: 

p = V 2 /47rGr 2 , (r<r h ); p = 0, (r > r h ), (56) 

where, as before, G is the gravitational constant. The velocity Vq is the (constant) circular velocity 
of the halo for r < r^, which we assume to be Vo = 220 kms" 1 as in Section There exist 
other parameterizations for the mass density that produce a smoother truncation at (e.g., Lin 
Ik Lynden-Bell 1982; Wilkinson & Evans 1999), but use of such a profile would not change the 
basic argument presented here. The gravitational potential is 

^ = V 2 [l+\n{r h /r)l (r <r h ); * = V 2 (r h /r), (r > r h ), (57) 

The enclosed mass is 



M(r) = rV 2 /G, (r<r h ); M(r) = r h V 2 /G = M MW , (r>r h ), 



(58) 



-31- 



where Mmw is the total mass of the Milky Way. It is useful to define some ancillary quantities: 

C = ^Kmc, i = \{v h Mc/V G f, M' = r LMC vl MC /2G. (59) 

The condition £ > that the LMC is bound to the Milky Way implies for the mass enclosed 
within the present radius of the LMC 

M(r L Mc)>M', (C<1); M(r LMC ) > M'/(l + In C), (C > !)• (60) 

Correspondingly, the total mass of the Milky Way must satisfy 

M MW >M', (C<1); M MW >M'C/(l + mC), (C > !)• (61) 

These equations illustrate a basic property of the models that fit the data: one can decrease the 
required mass M(rLMc) inside tlmc by increasing (, but this increases the total mass Mmw- Over 
the domain £ > 1, the function h(Q = £/(l + ln£) has a minimum h(Q = 1 at £ = 1. So if one 
has no knowledge of the extent of the halo, as quantified by then the strict lower limit on the 
mass of the Milky Way is 

M MW > M'. (62) 

In the present situation we do have a limit on £. We have assumed knowledge of Vb, and this fixes 
the normalization of the mass density (eq.(5(|). This allows one to calculate the minimum halo 
extent = (Ylmc for which Mmw = r hVo /G satisfies equation (|6l|). The result is: 

C>£, (£<i); C>^~\ (C>1). (63) 

Substitution of tlmc and vlmc from equations (54) and ( |55| ) and use of a simple Monte-Carlo 
scheme to propagate errors yields 

M MW > 4.3 x 10 11 M Q , r h > 39 kpc, (68.3% confidence). (64) 



It is interesting to compare these constraints from the assumption that the LMC is 
bound to the Milky Way with the results from studies of the ensemble dynamics of all 
the Milky Way satellites. Kochanek (1996) obtained for the mass of the Milky Way inside 
50 kpc that M(50 kpc) = (4.9 ± 1.1) x 10 11 M Q and Wilkinson & Evans (1999) obtained 
M(50 kpc) = 5.4^3g x 10 11 Mq. Both studies indicated that the dark halo extends considerably 
beyond 50 kpc and that Mmw exceeds M(50 kpc) by a factor of a few. These results are fully 
consistent with equation (]64]). 

The mass constraint in equation (|64|) is a factor ~ 2 smaller than the limit obtained by 
Pedreros et al. (2002), also from the LMC proper motion. Their result assumes that all the mass 
of the halo is enclosed within the present radius of the LMC. This is in contradiction with the 
results of Kochanek (1996) and Wilkinson & Evans (1999), and the results in equation ( |64| ) should 
therefore be more accurate. 
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9.5. Magellanic Stream 

The properties of the LMC orbit and the Milky Way dark halo are important parameters 
in models of the Magellanic Stream. Recent models generally assume an extended, usually 
isothermal, dark halo, consistent with the arguments of Section |9.4| . Models that fit the generic 
properties of the Stream suggest that the LMC is just past pericenter in its orbit. The predicted 
transverse velocity of the LMC has generally been found to be in the range vt = 400-500 kms -1 
approximately due East (reviewed by, e.g., Westerlund 1997). This is consistent with the observed 
transverse velocity (see Figure ||). In Galactocentric coordinates, all authors tend to agree that 
the component ULMC.rad is small. For fLMCtam a wide range of physical arguments has yielded 
values that include t>LMC,tan = 369 kms" 1 (Lin Sz Lynden-Bell 1982), 355 kms" 1 (Shuter 1992), 
352 kms" 1 (Heller & Rohlfs 1994), 339 kms" 1 (Murai & Fujimoto 1980), 320 kms" 1 (Liu 1992) 
and 285 kms" 1 (Gardiner et al. 1994; Gardiner & Noguchi 1996), respectively. The observationally 
derived value is t>LMC,tan = 281 ± 41 kms -1 (equation (|55|)), which is consistent with the lower 
range of the predicted values. The data are most consistent with the models of Gardiner et 
al. (1994) and Gardiner & Noguchi (1996), which are some of the most detailed models that have 
been constructed for the Magellanic Stream. This agreement is actually quite remarkable, given 
the many uncertainties and simplifications involved in models of the Magellanic Stream. The 
observational error on t>LMC,tan is almost small enough to start ruling out specific models. 

Lin et al. (1995) modeled the Magellanic Stream using a different approach than most 
authors. They fixed the orbit to that implied by the proper motion measurement of Jones et 
al. (1994), which was one of the first measurements to be made, and then explored different 
halo potentials. The low Galactocentric velocity t>LMC,tan = 213 ± 49 kms -1 implied by this 
proper motion measurement could only be reproduced with a rather 'lean' dark halo. With only 
~ 3 x 10 11 Mq inside 50 kpc their mass model is inconsistent with the results obtained and 



discussed in Section |9.4| . We attribute this to the limited accuracy of the Jones et al. (1994) proper 
motion measurement. The average proper motion inferred from all the data in Table Q is in fact 
consistent with canonical (isothermal) models of the dark halo. 

The Jones et al. value of VLMC.tan was consistent with the circular velocity of the Milky 
Way dark halo. Based on this it has often been said in the literature that the LMC is on an 
approximately circular orbit around the Milky Way. The results presented here do not support 
this view. They are consistent with the orbit advocated by Gardiner et al. (1994) and Gardiner &: 
Noguchi (1996), which has an apocenter to pericenter ratio of ~ 2.5 : 1. The perigalactic distance 
is ~ 45 kpc and the present orbital period around the Milky Way is ~ 1.5 Gyr. 

There remains considerable debate about the detailed origin of the Magellanic Stream. While 
tidal stripping effects have undoubtedly played an important role, it remains quite possible that gas 
dynamical effects and ram-pressure stripping have also played an important part in the shaping of 
the Stream (e.g., Heller & Rohlfs 1994; Moore & Davis 1994). In view of this it is unlikely that any 
modeling of the Stream will soon lead to unambiguous predictions for three-dimensional velocity 
of the LMC and its orbit. The only thing that can be safely concluded now is that observational 
measurements of the three-dimensional velocity of the LMC are generically consistent with models 
for the Magellanic Stream. 
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10. The Predicted LMC Proper Motion Velocity Field 

With future astrometric observations it may become possible to measure the proper motion 
velocity field of the LMC. This velocity field can be calculated with the formulae of Section ^ for 
any combination of model parameters. We have performed this calculation for our best model of 



the LMC structure and kinematics, characterized by the parameters given in equations ( |35| , 37 



40, 42, 48). The rotation velocity curve V(R') was parameterized as in equation (|3q), with the 
parameters 

V = 49.7kms -1 , Ro/D Q = 0.055, = 2.68, (65) 

chosen to best fit the points with error bars in Figure ||. In the absence of any observational 
information on dO/dt, this quantity was assumed to be zero. Figure |9] shows the resulting 
(residual) proper motion velocity field (i.e., the proper motion of the LMC CM was subtracted). 

The predicted proper motion velocity field is quite complex, and lacks any of the simple 
symmetries that one might have expected. While a small part of this is due to the fact that a plot 
of right ascension versus declination provides a distorted view of the sky, most of it is due to the 
fact that there are several different components that contribute to the velocity field. To provide 
some insight into the results, we show the individual components that contribute to the velocity 
field in Figure [n]. The top left panel shows the contribution due to the CM velocity component 
vcm (see eq. fL3|). This component dominates the total velocity field at most positions. It arises 
from the fact that one observes different components of the CM velocity vector at different 
positions on the sky. In particular, one observes part of the systemic velocity of the LMC in the V2 
component, so that v\ + / v\ unless one observes at the position of the CM (see eq. Jl3]]); here 
V2 and are the velocity components perpendicular to the line of sight, as defined in Section 0. 



The top right panel in Figure 10 shows the contribution due to the internal rotation component 
Vint (see eq. [21]) for the best model. This contribution has the expected morphology of a clockwise 
rotation in the plane of the sky. However, it does have a non-zero radial component (t>2 ^ 0) and 
a tangential component that varies as a function of position angle. The contribution of the di/dt 
component to the velocity field in Figure is small. However, it still is useful to have some insight 



into the morphology of this component. The bottom left panel of Figure 10 shows this component 



for an assumed value (di/dt) = 1 masyr -1 , which exceeds the value in our best fit model by a 
factor of ~ 3. For comparison, the bottom right panel of Figure 1C shows the contribution of a 
dQ/dt component of the same size (i.e., 1 masyr -1 ). As mentioned above, the latter contribution 
was set to zero in Figure ^ for lack of observational constraints. 

Figures |9]and 1C demonstrate a few important issues. First, as one moves away from the LMC 
center the residual proper motion velocity field can be quite significant as compared to the proper 
motion vector of the CM. Within the range of positions shown in Figures the predicted proper 
motion can differ by more than 1 masyr -1 from the value for the CM. Second, the majority of the 
residual proper motion is due to the velocity vector of the CM, and not due to the rotation of the 
LMC itself. It is therefore crucial that this contribution be properly included in any interpretation 
of proper motion measurements that are not made very close to the CM. And third, there may be 
contributions to the proper motion field due to precession and nutation of the LMC disk plane. 
These cannot be neglected in studies where high accuracy is desired. 
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11. Prospects for Kinematical Distance Determination of the LMC 

As described in Section [5^, there is still considerable uncertainty and controversy in the 
determination of the LMC distance Dq. Any new method that can shed additional information 
on this quantity is therefore potentially valuable. In this context it is interesting to note that 



equation (29) can be solved for the LMC distance Dq to obtain 

D = w ts /[fi s + (di/dt)}. (66) 

The quantity Wf s is determined by the line-of-sight velocity field, as is the line-of-nodes position 
angle G (eq.(37j). The proper motion component fj, s along the line of nodes can be measured. 
So if di/dt can be measured, or if one assumes that it is known, e.g., identical to zero, then 
this equation provides a new method for determining the LMC distance. The method is purely 
kinematical and bypasses all knowledge of stellar evolution, with its associated uncertainties. So 
it is complementary to most existing methods for determination of the LMC distance. 



The underlying idea of equation ( |66[ ) is that the line-of-sight velocity field yields the transverse 
velocity vector of the LMC in km/s, whereas proper motion measurements yield the same vector 
in milli-arcsec per year. Combination of the two measurements yields the distance. It is not 
a problem that the line-of-sight velocity field only constrains one component of the transverse 
velocity: one only needs one component to determine the distance. This method is essentially a 
variant to the 'radial velocity gradient' method which has been used to determine the distance 
to open clusters (e.g., Gunn et al. 1988). However, open clusters do not rotate appreciably, so 
there is no difficulty in disentangling the dynamics of the cluster from the apparent solid body 
rotation induced by the transverse motion. Gould (2000) was the first to point out that the 
method can also be used for the LMC, despite the fact that the LMC is itself rotating. Our 
approach is a variation on that proposed by Gould, but with important differences. First, we 
have devised a methodology and software that allow a practical application to real data. Second, 
Gould's suggested approach relied on the assumption that the kinematic line of nodes can be 
approximated to be the same as the photometric major axis. This is true for a circular disk, but 
is now known to be highly incorrect for the LMC (Paper II showed that the major axis position 
angle is PA ma j = 189.3° ± 1.4°, which differs from the kinematic line of nodes by ~ 60°). We 
have shown here that that the kinematic line of nodes can be determined independently from the 
line-of-sight velocity field, so that the method can still be used even if the photometric major axis 
does not trace the line of nodes. 

Two major hurdles must be overcome for the proposed method to ever become practical: 
first, accurate independent knowledge must be obtained on the value of di/dt; and second, the 
formal errors in our knowledge of wt s and fj, s must be reduced. These points are illustrated by 
Figure 11. It shows the 68.3% confidence region in the plane spanned by the parameters di/dt 
and m — M (the distance modulus corresponding to Dq), obtained with the help of equations ( f27| , 
37, 66). In Section [6| we used the independent estimate of the distance given by equation (f40|) 



to obtain an estimate for di/dt (eq. |4^]). If instead the goal is to determine m — M, then one 
needs an independent estimate of di/dt. This estimate must be accurate to 0.03 mas yr _1 to obtain 
m — M to better than 0.05 mag. Gould (2000) did not address the influence of di/dt, which is the 
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same as assuming that di/dt = 0. This is clearly not acceptable to obtain an accurate distance 
estimate. With the present data it implies that m — M = 19.08 ± 0.36, which is inconsistent with 
the generally accepted value (eq. |40[j ) at the 1.6a level. Even if di/dt were known with infinite 
accuracy, the method wouldn't presently be very competitive. At fixed di/dt, the error va. m — M 
is in the range 0.3-0.4 mag, too large to contribute meaningfully to the debate over the LMC 
distance. 

The main advance for application of the proposed method will come from future astrometric 
space missions such as SIM and GAIA. These can measure the proper motion of the LMC CM 
(and thus the component fi s that features in eq. pH]) to 1% or better (Gould 2000). This will 



reduce the width of the confidence band in Figure 11. In addition, they will yield information 
on the internal proper motion velocity fields of the LMC (for which predictions were shown in 
Figures |9] and [Tc| ) . The three components (v\, v%, V3) of the three-dimensional velocity (eq. 0]) 
each constrain a combination of the quantities Dq, di/dt and dQ/dt, as shown in Section These 
can be combined to obtain independent estimates of all three of these quantities. Information on 
the rotation of the LMC in the plane of the sky will also yield an estimate of the rotation curve 
V(R') in mas/yr, which can be combined with measurements obtained from line-of-sight velocities 
in km/s to obtain another constraint on the LMC distance. This is called 'rotational parallax' 
(Oiling & Peterson 2002). With sufficient information it may also become possible to constrain 
the influence of non-circular streamlines, which will reduce the effect of systematic errors on the 
results. 

To obtain an accurate distance measurement one will need not only improved proper motion 
data (to reduce the uncertaintity in fi s in eq. f66|| ), but also improved information on the 
line-of-sight velocities (to reduce the uncertaintity in wt s in eq. J66[). It should be possible to 
obtain velocities for ~ 10 times more carbon stars than are currently available, based on the total 
number of carbon stars in the LMC (Kontizas et al. 2001; Demers, Dallaire & Battinelli 2002; 
Alves et al., in prep.). This will reduce the errors in wt s and by a factor of ~ 3. The internal 
proper motion velocity fields that can be obtained with future astrometric missions can be used to 
further reduce the error in O. 



12. Summary and Conclusions 

Disk galaxies form an important constituent of the Universe and it is crucial to have an 
adequate understanding of their kinematics. The nearest disk galaxy is our own Milky Way. Our 
position in the Milky Way has allowed us to learn many things about it that are inaccessible for 
other galaxies. However, there are many things about the kinematics of our Milky Way that are 
understood less well than for disk galaxies in general (e.g., the behavior of the rotation curve 
outside the solar radius). This is due to the fact that we reside inside the Milky Way, so that 
the Milky Way spreads over the entire sky. This necessitates the study of Milky Way kinematics 
in terms of (l,b,v) diagrams, which are notoriously complicated to interpret (e.g., Binney & 
Merrifield 1998). By contrast, almost all other disk galaxies, including nearby ones such as 
M31 and M33, are so far away that they span <J 1° on the sky. This allows one to model their 
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kinematics fairly accurately under the assumption that the sky has no curvature over the area of 
the galaxy. The LMC is unique in that it is only disk galaxy that falls between these extremes: it 
is not close enough to cover the whole sky, yet close enough to have its observables influenced by 
the curvature of the sky. 

To fully understand the implications of the large extent of the LMC on the sky, we have 
derived general expressions for the velocity fields of the LMC, both along the line-of-sight and 
in the plane of the sky. We have analyzed these expressions to understand which quantities are 
uniquely constrained by the data, and which other quantities are degenerate. While we are not 
aware of any previous presentation in the literature of our general expressions for the LMC velocity 
field, this does not mean that previous authors have failed to see the importance of the large 
angular extent of the LMC. The most importance consequence is that the transverse motion of the 
LMC introduces a spurious solid-body rotation component in the observed line-of-sight velocity 
field, and this has certainly been realized by nearly all authors in the past decades. However, most 
authors have treated this component as a 'contamination', which needs to be subtracted before 
the internal kinematics of the LMC can be studied (e.g., Kim et al. 1998; Alves & Nelson 2000). 
While this is a valid viewpoint, it requires an accurate measurement of the LMC proper motion, 
which has not generally been available. In addition, this approach fails to realize that there is 
important information in the velocity field that can be extracted from the data without requiring 
prior knowledge of the proper motion. While this has been recognized previously by some authors, 
interpretations were always accompanied by the assumption that the line of nodes coincides with 
the photometric major axis of the LMC (e.g, Feitzinger et al. 1977; Meatheringham et al. 1988; 
Gould 2000). While this may have seemed reasonable (it is true for a circular disk), this is now 
known to be invalid (Papers I and II). The present paper represents the first exposition of what 
assumptions are actually necessary and/or reasonable, and what can actually be learned from the 
observed line-of-sight velocity field with and without these assumptions. 

We have fitted the general expression for the LMC line-of-sight velocity field to data for 1041 
carbon stars. Several things can be inferred even without any knowledge of the LMC transverse 
velocity and its associated solid-body rotation component. Most importantly, the position angle 
of the line of nodes is uniquely constrained at = 129.9° ± 6.0°. This is consistent with the 
value = 122.5 ± 8.3° determined geometrically in Paper I. The observed drift in the center of 
the isophotes at large radii is also consistent with these estimates, when interpreted as a result of 
viewing perspective (as discussed in Paper II). With three independent arguments there can now 
be little remaining doubt that the position angle of the line of nodes is in fact quite different from 
the major axis position angle, PA ma j = 189.3° ± 1.4°, with the corollary that the LMC must be 
intrinsically elongated (Paper II). 

The LMC dynamical center is determined by the data to ~ 0.4° accuracy in each coordinate. 
The inferred position is consistent with the center of the bar and the center of the outer isophotes. 
However, it is offset by 1.2° ± 0.6° from the kinematical center of the HI. The most plausible 
interpretation of this is that the structure and kinematics of the HI in the LMC are highly 
disturbed. This is quite consistent with our general knowledge of the gas morphology in the 
Magellanic Clouds system. Almost all previous studies of the kinematics of the LMC, even 
the ones using stellar tracers (carbon stars, planetary nebulae, etc.), have generally fixed the 
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kinematical center a priori to coincide with the HI kinematical center. The results of the present 
study indicate that this not appropriate. 

The LMC rotation curve can only be inferred with knowledge of vt c , the transverse velocity 
component along the line of nodes. This component is itself not (strongly) constrained by the 
velocity field. By contrast, the orthogonal component Vf S is constrained by the velocity field, but 
only through the quantity wt s = vt s + Do(di/dt). We find that wt s = —402.9 ± 13.0 kms -1 . This 
provides a non-trivial constraint on a combination of three apparently unrelated quantities: the 
transverse velocity component of the LMC in the direction perpendicular to the line of nodes, the 
LMC distance, and the rate of inclination change di/dt. 

To further interpret the results, we have compiled various LMC proper motion measurements 
from the literature. For two of the measurements, the ones by Jones et al. (1994) and Pedreros et 
al. (2002), we have calculated an improved correction for the offset of the field under study from 
the LMC CM. The weighted average proper motion, combined with an estimate for the LMC 
distance, yields both vt c and vt s . 

The vts inferred from the proper motion data can be combined with wt s inferred from the 
velocity field analysis and an estimate of the LMC distance to determine the rate of inclination 
change of the LMC disk: (di/dt) = —0.37 ± 0.22 masyr -1 = —103 ± 61 degrees/Gyr. This is not 
inconsistent with what can reasonably be explained as a result of precession and nutation of the 
LMC disk due to tidal torques from the Milky Way (Weinberg 2000). 

The vtc inferred from the proper motion data can be combined with the results from 
the velocity field analysis to yield the LMC rotation curve. We find that it flattens out at 
V = 49.8 ± 15.9 kms^ 1 . The error is dominated by our knowledge of the LMC proper motion. 
The rotation curve amplitude inferred here is ~ 40% lower than what has been indicated by all 
previous analyses of the LMC kinematics, including, e.g., the carbon star study of Alves & Nelson 
(2000) and the HI study of Kim et al. (1998). This is due to the fact that these studies used or 
imposed values of v tc , vt s , wt s , di/dt and/or that are inconsistent with the data and the analysis 
presented here. 

Upon correction for asymmetric drift, the circular velocity of the LMC is estimated to be 
^circ = 64.8 ± 15.9 kms" 1 . This implies an enclosed mass inside the last measured data point 
M LM c(8.9 kpc) = (8.7 ± 4.3) x 10 9 M . This corrects some previous estimates, which have ranged 
to values as high as 2 x 10 10 M . The inferred mass exceeds the sum of the estimated masses of 
the stellar disk (~ 2.7 x 10 9 M Q ) and the gaseous disk (~ 0.5 x 10 9 M ). This implies that the 
LMC is embedded in a dark halo, consistent with the relative flatness of the observed rotation 
curve out to the last measured data point. 

To estimate the tidal radius of the LMC we note that the traditional formulae for calculating 
this quantity are incorrect for galaxies that are embedded in a dark halo. We derive alternative 
formulae that are appropriate for the Milky Way - LMC system, which yield an LMC tidal radius 
r t = 15.0 ± 4.5 kpc (i.e., 17.1° ± 5.1°). As expected, this exceeds the radii out to which the LMC 
isopleths are observed to be regular. 

The carbon star analysis yields an average line-of-sight velocity dispersion a = 20.2±0.5km s^ 1 . 
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The radial profile of the velocity dispersion does not fall steeply enough with radius to be consistent 
with a constant scale-height disk. Application of the isothermal flared disk model formalism of 
Alves & Nelson (2000), in which the vertical density profile is proportional to sech 2 (,z/2!o)> implies 
that zq increases from 0.27kpc at the LMC center to 1.5 kpc at a radius of 5.5 kpc. More generally, 
the importance of rotational support and the disk scale height are both directly related to the 
quantity V/a. For the LMC carbon stars we find V/a ~ 2.9 ± 0.9. Compared to the Milky Way 
disk components, this is smaller than the values for both the thin disk (V/a ~ 9.8) and the thick 
disk (V/a ~ 3.9). Simple arguments for models stratified on spheroids indicate that the LMC 
disk could have an (out-of-plane) axial ratio of ~ 0.3 or larger. These results indicate that the 
LMC is a much thicker disk system than has been assumed previously. This is consistent with the 
predictions of iV-body simulations by Weinberg (2000) that address the Milky Way tidal influence 
on the LMC. 

To lowest order, the LMC self-lensing optical depth depends exclusively on the observed 
velocity dispersion and not separately on either the galaxy mass or scale height. Although we find 
a different rotation curve amplitude for the LMC than previous authors (which affects the mass 
and scale height of the LMC), we do find the same velocity dispersion. It is therefore not expected 
that the new insights from the present paper will substantially alter existing estimates for the 
LMC self-lensing optical depth. 

The proper motion data yield an estimate of the transverse velocity of the LMC. This velocity 
is similar to what can be estimated completely independently from the line-of-sight velocity field 
(for reasonable values of di / dt) . This provides confidence in the proper motion data, and indicates 
that there is no reason to believe that they are plagued by large systematic errors. The transverse 
velocity is vt = 406 km s -1 in the direction of position angle @t = 78.7°, with errors of ~ 40kms~ 1 
in each coordinate. This can be combined with the line-of-sight velocity inferred from the carbon 
star data, t> sys = 262.2 ± 3.4 kms -1 , to determine the LMC velocity vlmc i n the Galactocentric 
rest frame. After correction for the reflex motion of the sun we find that ulmc = 293 ± 39 kms -1 , 
with radial and tangential components t>LMC,rad = 84 ± 7kms~ 1 and t>LMC,tan = 281 ± 41 kins" 1 , 
respectively. This provides one of the most accurate space velocities of any Milky Way satellite. 
The result is at the low end of the range of velocities that has been predicted by models for 
the Magellanic Stream. The data are fully consistent with the models of Gardiner et al. (1994) 
and Gardiner & Noguchi (1996), which are some of the most detailed models that have been 
constructed. The implied orbit of the LMC has an apocenter to pericenter distance ratio ~ 2.5 : 1. 
The perigalactic distance is ~ 45 kpc and the present orbital period around the Milky Way is 
~ 1.5 Gyr. The constraint that the LMC is bound to the Milky Way provides a robust limit on the 
minimum mass and extent of the Milky Way dark halo. For an isothermal halo model truncated 
at r/i one finds that M MW > 4.3 x 10 11 M and r h > 39 kpc (68.3% confidence). 

In Figures [9| and 10 we have presented predictions for the LMC proper motion velocity field, 
based on our best fit model. This velocity field is dominated by the influence of the CM velocity 
(the proper-motion components of which vary with position in the LMC), and not by the internal 
rotation velocity component. The proper motions due to di/dt and dQ/dt play a minor, but 
non-negligible role. 
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The results of the line-of-sight velocity field analysis (in particular wt s and the position angle 
of the line of nodes) can be combined with proper motion data to obtain an estimate of the 
LMC distance. We argue that this method has promise, but show that it will not be competitive 
until: (a) di/dt can be accurately determined from independent data (e.g., measurement of the 
LMC proper motion velocity field); (b) the proper motion of the LMC CM is determined more 
accurately; and (c) a larger set of stellar line-of-sight velocities is obtained for analysis. 

The next advance in our understanding of LMC structure and kinematics will have to come 
from more accurate proper motion data. Such measurements will be possible with future space 
missions (e.g., SIM, GAIA), and can possibly be achieved with other means as well (e.g., a study 
by Alcock et al. is underway with the Advanced Camera for Surveys on Hubble Space Telescope). 
The uncertainties in our knowledge of many of the LMC's structural parameters are dominated 
by the uncertainties in the presently available proper motion measurements (Section [Tl]). This is 
true, e.g., for the LMC rotation curve V(R'), for its mass and tidal radius, and for di/dt. Our 
knowledge of these quantities will improve significantly with a more precise measurement of the 
proper motion of the LMC CM. This will also determine the space motion of the LMC to high 
accuracy. Combined with models for the Magellanic Stream this should allow the extraction of 
improved constraints on the density distribution of the Milky Way dark halo. Measurements of 
the proper motion velocity field of the LMC will yield better constraints on the LMC distance, 
by providing independent constraints on di/dt, dQ/dt, and possible deviations from circular 
streamlines. In addition, it will allow application of the method of rotational parallax, which 
compares the rotation amplitude in masyr -1 inferred from the proper motion velocity field with 
the same quantity in km/s inferred from the line-of-sight velocity field. 

We are grateful to Andrew Drake, Mario Pedreros and Mir Abbas Jalali for discussions and/or 
sharing results prior to publication. We dedicate this paper to the memory of Robert Schommer. 
He made many important contributions to the understanding of the structure and kinematics 
of the LMC, and he was a valued collaborator on the acquisition of carbon star velocity data 
analyzed in the present paper. 
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A. The Proper Motion of the LMC Center of Mass from Off-Center Measurements 

Some of the proper motion measurements that have been reported in the literature for the 
LMC pertain to fields that are at a considerable distance from the LMC CM. This includes in 
particular the measurements of Jones et al. (1994) and Pedreros et al. (2002). These measurements 
must be corrected for the orientation and rotation of the LMC disk to obtain an estimate of 
the proper motion of the LMC CM. These corrections have usually been made with the aid of 
kinematical models for the LMC that are less sophisticated than those presented here. Among 
other things, the position angle of the line of nodes that has generally been assumed is quite 
different from the value derived here and in Paper I. In view of this, we rederived the corrections 
that must be made to the measurements of Jones et al. (1994) and Pedreros et al. (2002) to obtain 
estimates of the proper motion of the LMC CM. The results are more accurate than the estimates 
quoted in these respective papers. 

The field studied by Jones et al. (1994) has J2000 coordinates a = 97.491° and 8 = -64.166°. 
The measured proper motion was pw = —1-26 ± 0.28 mas/yr and pjy = 0.26 ± 0.27 mas/yr. Given 
the CM position derived in Section | (eq. f37|), the field center has p = 8.298° and $ = 54.195°. 
With the inclination and line of nodes position angle given by equations (|35|) and (|37|), one obtains 
from equation (pi) that R'/D = 0.16. At this position in the LMC disk V(R') « 50 kms" 1 , 



according to Figure |y. We take v sys from equation ( p7| ) and Dq from equation ( |40| ) . The formulae 
of Section || then allow calculation of the proper motions (pw,^n) at the Jones et al. field center, 
for any assumed proper motion (pw,o, I^Nfi) of the LMC CM. At any value of (/iw,OiMJV,o) we 
choose di/dt so as to reproduce the value wt s that was inferred from the line-of-sight velocity field 
(eq. [37|). We determined those values of {pwfl-, MiV.o) that reproduce the Jones et al. measurement, 



which yielded pw,o = —1-36 ± 0.28 mas/yr and p^fi = —0.16 ± 0.27 mas/yr. 

The field studied by Pedreros et al. (2002) has J2000 coordinates a = 74.989° and 
5 = —64.378°. The measured proper motion was pw — —1-8 =t 0.2 mas/yr and pm = 0.3 ± 0.2 
mas/yr. We proceeded similarly as for the Jones et al. data. The field center has p = 6.111° and 
$ = 330.698°, and R'/D = 0.11. At this position in the LMC disk V(R') « 50kms _1 , according 
to Figure ©. The proper motion of the LMC CM implied by the Pedreros et al. measurement is 
^w,o = —1-83 ± 0.20 mas/yr and pn,o = 0.66 ± 0.20 mas/yr. 

The sizes of the corrections that must be applied to the Jones et al. (1994) and Pedreros 
et al. (2002) data to obtain an estimate of the proper motion of the LMC CM are 0.43 mas/yr 
and 0.36 mas/yr, respectively. At the distance of the LMC this corresponds to 103 kms -1 and 
86kms~ 1 , respectively, which considerably exceeds the rotation amplitude of the LMC. The 
reason for this is that the residual proper motion field (i.e., the proper motion field with the 
proper motion of the CM subtracted) is dominated by the velocity component of the LMC CM, 



and not by the velocity component due to internal rotation. This was discussed in Section 10. 
Because only a small part of the corrections is due to the rotation of the LMC, the corrections are 
quite insensitive to possible errors in the parameters that characterize the model. As a result, the 
formal errors in {pwfl-, ^Nfl) are dominated by the formal errors in the observed (pw^n)- 



The corrections that were derived assume that the position angle of the line of nodes is not 
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time dependent. If instead d@/dt ^ 0, then the v% term in equation (|i~6| ) has an extra component. 
However, the estimated (nwfi, I^Nfi) do not depend sensitively on dQ/dt. If dQ/dt = 1 masyr -1 , 
then (fJ,w,o, A'iV.o) changes by (—0.076,0.109) masyr -1 for the Jones et al. measurement and by 
(—0.093, —0.053) for the Pedreros et al. measurement. These corrections are small compared to 
the formal errors in the data. The iV-body simulations of Weinberg (2000) show that it is unlikely 
that the size of dQ/dt is much larger than 1 masyr -1 = 277.78°/ Gyr. 
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Fig. 1. — Illustration of the projected view of the sky. All vectors and angles lie in the plane of 
the paper (perspective was not used in the drawing). All quantities are rigorously defined in the 
text. The galaxy center of mass (CM) with coordinates (ao,5o) is chosen to be the origin O of 
an (x, y, z) coordinate system. The z-axis (not drawn) points vertically out of the paper, towards 
the observer. The angles p and = cf) — 90° define the projected position on the sky of a tracer 
with coordinates (a, S). The angle T is the angle between the (proper motion) vectors nw and 
towards the local directions of West and North, and the velocity components V2 and defined in 
the text. 
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OBSERVER' 




Fig. 2. — Schematic 'side view' of the observer-galaxy system. All quantities are rigorously defined 
in the text. The galaxy center of mass (CM) is the origin of an (x, y, z) coordinate system. The 
z-axis points towards the observer. In reality, the x-axis points vertically out of the paper but 
is shown here using drawing perspective. The distance from the observer to the CM is D$\ the 
distance from the observer to a tracer is D. The angle from the CM to the tracer, as seen by the 
observer, is p. The velocity of the tracer can be decomposed into three orthogonal components, 
vi, t>2 and U3, as described in the text. These components are shown schematically using drawing 
perspective. The component v\ = v\ os lies along the line of sight, and points away from the observer. 
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(South) 

Fig. 3. — Illustration of the observer's view of the galaxy disk. All quantities are rigorously defined 
in the text. The cartesian (x, y, z) coordinate system is the same as in Figure |], with the z-axis 
(not drawn) pointing vertically out of the paper, towards the observer. The (x',y',z') system is a 
second cartesian coordinate system, the axes of which are shown with drawing perspective. The 
(a/, y') plane is the plane of the galaxy disk, and the z'-axis is its symmetry axis. The plane of the 
galaxy disk is titled diagonally out of the paper. The inclination % is the angle between the (x, y) 
plane of the sky, and the (x',y') plane of the galaxy disk. The x'-axis is defined to lie along the 
lines of nodes (the intersection of the (x,y) plane of the sky, and the (x',y') plane of the galaxy 
disk). The angle = — 90° is the position angle of the line of nodes in the plane of the sky. 
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Fig. 4. — Illustration of the quantities used to describe the transverse velocity of the LMC CM. 
The transverse velocity vector has size vt and is directed along position angle @t = Of — 90°. It has 
components v x and v y towards the West and North, respectively. For analysis of the velocity field 
of the LMC it useful to decompose the transverse velocity vector into a sum of two orthoghonal 
vectors, as illustrated. One vector lies along the position angle of the line of nodes, and has 
length vtc- The other lies along the position angle — 90°, and has length —vt s (i-e., it has length 
vts in the direction + 90°). The velocity components shown in the figure correspond to the proper 
motions listed in parenthesis (upon division by the LMC distance Dq). 
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Fig. 5. — Carbon star line-of-sight velocity data from Kunkel et al. (1997) and Hardy et al. (2002, in 
prep.), as function of position angle on the sky. The displayed range of the angle $ is 0°-720°, so 
each star is plotted twice. Each panel corresponds to a different range of angular distances p from 
the LMC center, as indicated. The curves show the predictions of the best-fit model (calculated at 
the center of the radial range for the given panel). 
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Fig. 6. — (Top Panel) Rotation velocity V in the plane of the LMC disk as function of the 
cylindrical radius R' (expressed in units of Dq, the LMC distance; Dq ~ 50 kpc). The solid curve 
corresponds to the best fit to the data (see Figure ||) for a model parameterized by equation (p6| ) . 
However, each of the dashed curves provides the same acceptable fit, but for a different assumed 
value of the component vt c of the transverse velocity along the line of nodes. The values of vt c (in 
km/s) for the different curves are labeled. The radial range that is displayed corresponds to the 
radial range p < 13° for which kinematical data were included in the fit. The exponential disk scale 
length of the LMC is R' d « 1.4 kpc (Paper II), which corresponds to R' d /D m 0.028. The data 
points with error bars show the results of model fits to individual rings in the LMC disk plane, as 
described in Section fL% . These fits had V and as only free parameters; the quantity vtc was fixed 
to the value 253 km s -1 indicated by proper motion data (eq. (^3|), and the other model parameters 
were fixed to the values (eq. |p7|1 ) determined for the best-fit model shown in Figure || (Middle 
Panel) Radial profile of the line-of-sight velocity dispersion. (Bottom Panel) Radial profile of 
the line-of- nodes position angle G. 
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Fig. 7. — The (hw^n) plane spanned by the proper motion components of the LMC center of 
mass. The data points with error bars show available proper motion measurements, as listed in 
Tabel |]. Dotted ellipses are the corresponding 68.3% confidence regions (note that the l-a error 
bars for a two-dimensional Gaussian only contain 39.3% of the cumulative probability). The heavy 
solid ellipse is the 68.3% confidence region that corresponds to the weighted average of the data. 
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Fig. 8. — The (v x ,v y ) plane spanned by the transverse velocity components of the LMC center 
of mass. The heavy solid ellipse is the 68.3% confidence region obtained from the proper motion 
data listed in Table [j] and shown in Figure [7|. The dotted trapezoid is the 68.3% confidence region 
obtained from the LMC line-of-sight velocity field under the assumption that di/dt = 0. Along 
the long axis of the trapezoid all points are equally probable; along the short axis the probability 
distribution is Gaussian. The solid trapezoid shows the corresponding result if di/dt is chosen 
to best agree with the proper motion data. The arrows indicate the directions along which the 
components vt c and vt s of the transverse velocity are measured. The former is measured along 
the line of nodes, and the latter is measured perpendicular to the line of nodes. A change in the 
assumed di/dt corresponds to a shift of the trapezoid along the vt s direction. 
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Predicted LMC Proper Motion Velocity Field 




Fig. 9. — The proper motion velocity field predicted by our best model for the LMC structure and 
kinematics. At each position on in a grid in right ascension and declination we have calculated the 
predicted proper motion components {hw-,^n) using the formulae of Section |2[ The proper motion 
vector is shown as a line segment, which starts at the grid point for which it was calculated (thin 
solid dots). The proper motion of the LMC CM, which is shown as a heavy line segment starting 
at the LMC CM (heavy dot), was subtracted before plotting (in other words, the actual predicted 
proper motion at a given point is the vector sum of the line segment shown at that point and the 
heavy line segment). Each segment is shown with the nw component horizontally to the right and 
the fj,jsr component vertically upwards, with equal scales in mas/yr in both directions. The LMC 
CM proper motion vector in the model has length 1.71 masyr -1 , which determines the scale of all 
line segments. 
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Predicted LMC Proper Motion Velocity Field 

Precession/Nutation Component Only (di/dt=0 ; d0/dt=l mas/yr) 
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Fig. 10. — Vector components of the predicted LMC proper motion velocity field, displayed as in 
Figure (Top Left) The proper motions due to the CM velocity component vqm (see eq. [ 13 1 ) 
for the model shown in Figure ^. (Top Right) The proper motions due to the internal rotation 
component v- mt (see eq. fl2l|| ) for the model shown in Figure ||. (Bottom Left) The proper motions 
due to the disk precession/nutation velocity component v pn (see eq. |l6|]) for (di/dt) = 1 masyr -1 
and (dO/dt) = 0. (Bottom Right) The proper motions due to the disk precession/nutation 
velocity component -Upn for (di/dt) = and (dQ/dt) = 1 masyr -1 . For comparison, the model 
shown in Figure | has (di/dt) = -0.37 masyr 1 and (dQ/dt) = 0. 
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Fig. 11. — The distance modulus of the LMC m — M (in magnitudes) versus the rate of inclination 
change di/dt. The analysis of the line-of-sight velocity field and the available proper motion data 
constrain a combination of these parameters. The heavy solid curve shows the best fit, while the 
thinner surrounding curves indicate the la confidence region. The dashed horizontal line indicates 
the canonical distance modulus of the LMC inferred from other studies: m — M = 18.50 (see 
Section 5.2); the uncertainty on that distance modulus is ~ 0.1 mag. 
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Table 1. Measurements of the Proper Motion of the LMC Center of Mass 



paper 


label 






comments 






( masyr -1 ) 


( mas yr _1 ) 




Kroupa et al. (1994) 


PPM 


-1.30 ±0.60 


1.10 ±0.70 




Jones et al. (1994) 


JKL 


-1.36 ±0.28 


-0.16 ±0.27 


see Appendix [A| 


Kroupa & Bastian (1997) 


HIP 


-1.94 ±0.29 


-0.14 ±0.36 




Pedreros et al. (2002) 


P02 


-1.83 ±0.20 


0.66 ±0.20 


see Appendix [A| 


Drake et al. (2002) 


MACHO 


-1.40 ±0.40 


0.38 ±0.25 





Note. — Measurements of the proper motion of the LMC CM from various sources, as listed in 
column (1). The label used to indicate the measurements in Figure [7] is listed in column (2). The 
proper motions listed in columns (3) and (4) are characterized by the components in the directions 
towards the West {^w = — n a cos5) and North (/Utv = /i^), respectively. Listed uncertainties are 
formal la measurement errors. Jones et al. (1994) and Pedreros et al. (2002) reported proper 
motion data for fields at many degrees from the LMC CM. These measurements were used to 
estimate the proper motion of the LMC CM as discussed in the Appendix. These estimates are 
more accurate than those quoted in the original papers. 
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Table 2. Internal Kinematics of the LMC 



Bf/Do 


V 


Ay 


a 


Act 





A6 




( kms -1 ) 


( kms -1 ) 


( kms -1 ) 


( kms -1 ) 


(deg) 


(deg) 


0.009 


-27.9 


4.5 


14.2 


1.3 


96.4 


11.8 


0.028 


14.2 


4.1 


21.5 


1.4 


151.6 


37.7 


0.044 


25.2 


4.0 


19.7 


1.1 


120.9 


21.8 


0.060 


35.7 


5.9 


22.0 


1.4 


140.8 


12.0 


0.081 


57.3 


4.2 


17.0 


1.1 


141.2 


7.9 


0.096 


50.0 


3.1 


16.3 


0.9 


110.5 


7.0 


0.113 


39.4 


5.3 


16.3 


0.9 


142.4 


14.0 


0.130 


46.6 


8.0 


16.3 


1.4 


108.7 


9.5 


0.148 


55.8 


11.6 


22.1 


2.3 


132.0 


12.1 


0.178 


32.1 


16.5 


21.3 


2.2 


123.2 


38.7 



Note. - Radial profiles for the LMC of the rotation velocity V(R'), line-of-sight velocity 
dispersion a(R'), and kinematical line-of-nodes position angle Q(R'), as calculated in Section ^ 
and shown in Figure |6|. All quantities are listed with formal la random errors. The radius R' is 
the cylindrical radius in the disk, which in column (1) is expressed in units of the LMC distance 
Dq (which equals 50.1 kpc for m — M = 18.50). The quantity V is the in-plane velocity as defined 
in Section ^. 



